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ABSTRACT 

We present a sample of 723 optically selected BL Lac candidates from the SDSS DR7 spectroscopic 
database encompassing 8250 deg 2 of sky; our sample constitutes one of the largest uniform BL Lac 
samples yet derived. Each BL Lac candidate has a high-quality SDSS spectrum from which we 
determine spectroscopic redshifts for ^60% of the objects. Redshift lower limits are estimated for the 
remaining objects utilizing the lack of host galaxy flux contamination in their optical spectra; we find 
that objects lacking spectroscopic redshifts are likely at systematically higher redshifts. Approximately 
80% of our BL Lac candidates match to a radio source in FIRST/NVSS, and ~40% match to a 
ROSAT X-ray source. The homogeneous multiwavelength coverage allows subdivision of the sample 
into 637 radio-loud BL Lac candidates and 86 weak-featured radio-quiet objects. The radio-loud 
objects broadly support the standard paradigm unifying BL Lac objects with beamed radio galaxies. 
We propose that the majority of the radio-quiet objects may be lower- redshift [z < 2.2) analogs to 
high-redshift weak line quasars (i.e., AGN with unusually anemic broad emission line regions). These 
would constitute the largest sample of such objects, being of similar size and complementary in redshift 
to the samples of high-redshift weak line quasars previously discovered by the SDSS. However, some 
fraction of the weak-featured radio-quiet objects may instead populate a rare and extreme radio- weak 
tail of the much larger radio-loud BL Lac population. Serendipitous discoveries of unusual white 
dwarfs, high-redshift weak line quasars, and broad absorption line quasars with extreme continuum 
dropoffs blueward of rest-frame 2800 A are also briefly described. 

Subject headings: BL Lacertae objects:general — galaxies: active — quasars: general — surveys 



1. INTRODUCTION 

BL Lacertae objects compose an especially rare sub- 
class of Active Galactic Nuclei (AGNs) traditionally in- 
terpreted as low-luminosity radio galaxies with a rel- 
ativistic jet pointed toward the observer (e .g., see 
IBlandford fc Reei 11971 fUrrv & Padovanilll995D . Emis- 
sion from the relativistically boosted jet dominates the 
observed flux, yielding featureless (or nearly feature- 
less) optical spectra, strong radio and X-ray emission, 
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strong polarization, flat radio spectra, and rapid variabil- 
ity across the entire electromagnetic spectrum (e.g., see 
lKollgaardl[l99l IPerlman et al.l l200lh . Given their odd 
spectral characteristics, and the need to obtain optical 
spectra to confirm a candidate as a BL Lac object, the 
most efficient BL Lac selection algorithms tend to search 
for optical count erparts to strong r adio and/or X-ray 



sources (e.g., see IStickel et al.l Il991t IStocke et alj 
IPerlman et al.l 119981: ILaurent-Muehleisen et al.l 
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1999 



Landt et al.1120011: lAnderson et al.ll2007t IPadovani et all 
20071: iPiranomonte et all 120071 : iTurriziani et all 120071 : 
Plotkin et al.ll2008ft . 

Optically selected BL Lac samples have historically 
been extraordinarily difficult to assemble, but they are 
highly desired. Systematic searches based only on opti- 
cal characteristics would allow the exploration of BL Lac 
properties with minimal concern for biases introduced 
by radio and X-ray surveys (which generally tend to 
have shallower flux sensitivities than optical surveys.) 
Also, optically selected samples might reveal new pop- 
ulations of BL Lac objects, and they might even lead to 
serendipitous discoveries of other, perhaps new, classes 
of featureless objects. Unfortunately, the odd spectral 
characteristics of BL Lac objects conspire to make se- 
lection based solely on optical colors ineffective. For 
example, only 4 BL Lac objects were identified in 
the P alomar- Green Sur vey for ultraviolet-excess objects 
(PG. iGreen et all IT986) : another 3 PG sources initially 
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misclassified as DC (i.e., featureless) white dwarfs were 
identified as BL La c objects only after the ROSAT All 
Sky Survey (RASS.IVoges et allll99fll2000h came online 
([Fleming et al.lll993f ). 

Attempts have been made to exploit other distinctive 
BL Lac properties, like polarization and flux variability, 
to derive optically selected samples. Such techniques, 
however, require observing large areas of the sky at mul- 
tiple epochs just to assemble manageable lists of candi- 
dates suitable for spectroscopic follow-up. For example, 
it would seem that the strong polarization of BL Lac ob- 
jects (P > 2%) should allow for efficient identification of 
BL Lac candidates. However, mos t attempts along these 
lines have so far been unsuccessful (llmpey fc Brandlll982l ; 
iBorra fc Corriveaulll984t Uannuzi et al.lll993t ). probably 
because the polarized flux is variable with a low duty 
cycle: for example, X-ray selected BL Lac objects from 
the Einstein O bservatory Extende d Medium-Sensitivity 
Survey fEMSS. IStocke et al.lll991f) spen d over half their 
time with P < 4% (jJannuzi et al.lll99l ). 

Not until we entered the age of large-scale spectro- 
scopic digital sky surveys has optical selection of BL Lac 
objects been a realistic endeavor. Traditionally, the pri- 
mary limitation of efficient BL Lac recovery algorithms 
is the large number of false positives returned, even when 
radio and/or X-ray information is consulted. Optical 
spectroscopy is then required to make firm BL Lac iden- 
tifications. However, with a digi tal survey such a s the 
Sloan Digital Sky Survey (SDSS, ^rFeTaDEOOO), this 
task is less daunting because the requisite spectra al- 
ready exist. Furthermore, the most relevant spectra can 
be searched within databases, allowing for automated re- 
moval of many of the contaminants that previously had 
to be culled manually (after a significant investment in 
telescope time.) 

The first successful attempt at assembli ng an optically 
selected sam ple is the 2QZ BL Lac survey (Londish et al. 
[20021 [2007t ). derived from the 2-degree field (2dF) and 
the 6-degree field (6dF) quasi-stellar object Redshift Sur- 
veys (|Croom et al.ll2004l) . They searched over ~10 3 deg 2 
and recovered 7 confident BL Lac objects. All 7 ob- 
jects were additionally identified with radio and X-ray 
sources (post-selection), and they also showed optical 
flux variations in a photometric monitoring campaign 
over 2002-2004 (jNesci et al.ll2005[ ). An additional object 
(z = 0.49 4) is identified as a radio-quiet /weak BL Lac 
candidate (jLondish et al.|[2004f ). This source lacks both a 
radio and an X-ray detection (with a radio to optical flux 
ratio limit placing it firmly in the radio-quiet regime 12 ), 
but its optical properties are similar to (radio-loud) 
BL Lac objects. The recovery of only a single extragalac- 
tic object with a featureless spectrum lac king radi o emis - 
sion by the 2QZ survey confirms the lStocke et al.l (|1990l ) 
result that radio-quiet BL Lac objec ts must be extremel y 
rare if they even exist at all (also see Jann uzi et al.lll993f) . 

Searching 2860 deg 2 , iCollinge et all (|2005l ). hereafter 

12 Radio-quiet quasars are commonly defined as having radio 
to optical flux ra tios R < 10 (e.g., see Kcllcrmann et al,lll9891 ; 
Stockc ct al. 1992i). Throughout this paper, we adopt the approx- 
imately similar a ro < 0.2, where a ro is the radio-to-optical broad- 
band spectral index: a ro = — log(L o /L r )/5.08, where L and L r 
are optical and radio sp ecific luminosities at rest-frames 5000 A 
and 5 GHz, respectively (Stockc ct al. 1983). 



IC051 assembled a larger optically selected sample from 
the SDSS. Their sample contains 386 nearly feature- 
less spectra, including 240 probable BL Lac candidates 
and 146 possible BL Lac candidates (the latter are 
most likely DC white dwarfs, proper motions are con- 
sulted for their classification.) The vast majority of the 
probable candidates either have firm extragalactic red- 
shifts, or they match to a radio source in the Faint Im- 
ages of the Radio Sky at Twenty-cm survey (FIRST, 
iBecker et al1fl995l) and/or the NRAO VLA Sky Survey 
(NVSS. lCondon et al.l fl998). or to a RASS X-ray source 
(correlations to multiwavelength catalogs are performed 
p ost-se lection.) 

IC05I also present a list of 27 intriguing radio- 
quiet /weak BL Lac candidates. These objects have 
nearly featureless optical spectra, and their radio fluxes 
(or limits in the cases of radio non-detections) place 
their radio to optical flux ratios within or tantalizingly 
close to the radio-quiet regime. All but one of these 
also lack X-ray detections in RASS (but their extant X- 
ray limits are not sensitive enough to declare them X- 
ray weak.) Some of these objects are at high-redshift 
(z > 2.2), and those high-redshift objects may be alter- 
nately described as members of a population of high- 
redshift weak line quasars (WLQs) discovered b y the 



SPSS fe.g.. see iFan et al l Il999t [Anderson et al 



iShemmer et aLll2006Ll20'09TlDiamond-Stanic et al.ll2009f) 



2001; 



Three of these high-redshift objects w ere detected in 
follow-on Chandra X-ray observations bv IShemmer et ahl 
(2001), and they are X-ray weaker than typical radio- 
loud BL Lac objects. Other radio-quiet BL Lac candi- 
dates in IC05I have lower redshifts, and it is unclear if 
they are more closely related to BL Lac o bjects or to 
WLQs. Regardless of their true nature. IC05l demonstrate 
that optical selection of BL Lac objects from the SDSS 
is efficient, producing a large enough sample to reveal es- 
pecially rare sub-populations of objects with featureless 
spectra in numbers comparable to entire venerable radio 
and X-ray selected BL Lac samples. 

Here, we expand on IC05I and continue the search for 
BL Lac objects in the SDSS (which now covers almost 
three times the sky area as the IC05I sample) . We re- 
cover 723 objects. Our automated selection algorithm 
is described in SJ21 We visually inspect ^23,000 spectra, 
and we discuss the removal of various classes of contam- 
inants (including stars, post-starburst galaxies, quasars, 
etc.) in Sj3] we also present lists of serendipitously discov- 
ered unusual white dwarfs, WLQs, and unusual broad 
absorption line quasars. In |2] we describe the removal 
of contaminants that survive visual inspection (including 
featureless white dwarfs and large red galaxies.) The fi- 
nal sample is presented in $5j and we describe a spectral 
decomposition of the BL Lac host galaxy from the AGN. 
We use this decomposition to measure optical spectral 
indices for the AGN, to estimate lower redshift limits by 
assuming BL Lac host galaxies are standard candles, and 
to calculate decomposed AGN component optical fluxes 
and luminosities. We also correlate the 723 BL Lac can- 
didates to the FIRST and NVSS radio surveys, and to 
the RASS X-ray survey. The sample is discussed in fj6l 
which focuses on a particularly interesting subset of 86 
radio-quiet objects with weak spectral features. Finally, 
our main conclusions are summarized in [J71 Throughout, 
we use quasar to refer to any AGN (regardless of lumi- 
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nosity or radio loudness and including BL Lac objects), 
and we adopt H = 71 km s" 1 Mpc -1 , fi m = 0.27, and 
Oa = 0.73. 

2. SAMPLE SELECTION 

The optically selected BL Lac sample is taken from 
8250 deg 2 of the S PSS Data Release 7.1 (DR7.1, 
lAbazaiian et alj I2009D spectroscopic survey, which in- 
cludes 1.4 x 10 6 spectra. The SDSS is a multi- 
institutional effort to image 10 4 deg 2 of the north Galac- 
tic cap in 5 optical filters covering 3800 to 1 0,000 A (e.g., 
see lFukugita et al.lll996t IGunn et al.lll998l) . with follow- 
up moderate resolution spectroscopy (A/AA ~ 1800) of 
10 6 galaxies, 1 5 quasars, and 10 5 unusual stars (e.g., 
sec lYork et al.ll2000h . Data are taken with a special pur- 
pose 2.5-mete r telescope locate d at Apache Point Ob- 
servatory (see IGunn et al. 2006), with astrometric accu- 
racy at the ^100 milli - arcse c level at the survey limit 
of r ~ 22 (jPier et al.1 12003D and typical photometric 
precision a t brighter magnit udes of approximately 0.02- 
0.03 mag (Ive zic et al.l 120041 ); 640 simultaneous spectra 
are obtained over a 7 deg 2 field with a multi-fiber optical 
spectrograph. Further technical details can be found in 
IStoughton etai] (12002) . The SDSS DR5 Quasar Catalog 
(|Schneider et al.ll2007f ) contains 77,429 quasars with reli- 
able spectroscopic redshifts. Our selection process iden- 
tifies a very small fraction of objects (<0.1%) in the DR5 
Quasar Catalog as BL Lac objects. Although the DR5 
Quasar Catalog does not formally include BL Lac ob- 
jects, some overlap is expected between our BL Lac sam- 
ple and the quasar catalog, as our BL Lac definition can 
include some very weak emission line objects. 

Recent SDSS data releases contain an increasing num- 
ber of stellar spectra. Around 10% of the spectra in 
DR7.1 were targeted by the Sloan Extensio n for Galactic 
Unde rstanding and Exploration (SEGUE, lYannv et al.l 
I2009D survey. SEGUE complements the SDSS Legacy 
survey (which focused on galaxies and quasars) , and it is 
designed to acquire spectra of 240,000 stars in the Galac- 
tic disk and spheroid. 

SDSS data releases typically contain improvements to 
the reduction pipelines. Many of the recent spectroscopic 
improvements pertain to better wavelength calibration 
(e.g., better flat fields, more robust algorithms to cali- 
brate wavelengths from arc lamps, etc.), as required to 
meet the SEGUE science goals. Wavelength calibration 
improvements are typically on the order of ~10 km s _1 
and should not impact BL Lac selection. However, an 
important improvement to the spectrophotometric cal- 
ibrati on was implemented starting wi th Data Release 6 
(DR6. lAdehnan-McCarthv et al.|[2007h . The flux scale of 
SDSS spectra is now tied to point spread function (psf) 
magnitudes of standard stars, instead of fiber magnitudes 
as in Data Releases 1-5. The effect is generally small over 
the wavelength coverage of the SDSS spectrograph, but 
flux densities in the far blue end of SDSS spectra can 
differ by as much as 30% before and after D R6 (e.g., see 
Figure 5 of lAdelman-McCarthv et al.H2007l ). 

The new flux scale can change our measures of blue 
spectral features, especially the strength of the Ca II H/K 
break (C, see Equation [1] below) . The BL Lac sample 
presented here is the first SDSS sample selected post- 
DR6, so we do not expect all previously identified SDSS 



BL Lac objects to make it into this sample (even if 
we were to apply identical selection approaches.) How- 
ever, the majority of previously identified SDSS BL Lac 
objects should still be recovered. Previously identified 
SDSS BL Lac objects that are now rejected because of 
the improved spectrophotometric calibrations should pri- 
marily be those with relatively large H/K breaks (i.e., 
0.3 < C < 0.4) as measured in their pre-DR6 spectra 
(but with C > 0.4 in their DR7.1 spectra.) We describe 
the strength of the Ca II H/K break, C, as the fractional 
change in continuum blueward and redward of rest-frame 
4000 A (e.g., see lLandt et"aT1l2002l ): 

n _ (U,r) - (U,b) „ 1 . nQR ( (h,r) ~ (h,b) \ 

C (U) = ' 14 + ' 86 l (h,r) ) ' 

(1) 

where (f„,b) and {f v , r ) refer to the average fluxes per 
unit frequency just blueward (3750-3950 A) and just 
redward (4050-4250 A) of the H/K break, respectively. 
The quantities (f\ t b) and (fx,r) similarly refer to the 
average fluxes per unit waveleng th. A typical ellipti- 
cal galaxy has C ~ 0.5 (e.g., see iDressler fc Shectmanl 
l!987tlMarcha et "aL1ll996f ) ; in the context of the standard 
BL Lac unification model, more highly beamed BL Lac 
objects h ave weaker H/K b reaks (i.e., smaller C mea- 
sures, see lLandt et al.ll200l . 

2.1. Selection Criteria 

We generally require each object's SDSS spectrum to 
pass the following criteria for inclusion in this BL Lac 
sample: (1) emission lines show rest-frame equivalent 
widths (REW) < 5 A; and (2) the Ca II H/K depres- 
sion has measured C < 0.4. These criteria are commonly 
adopted in the literature and have similarly been a pplied 
to assemble previo us SDSS BL Lac samples (e.g., \CM 
iPlotkin et aLll2008D . 

In general, we do not explicitly require a close match 
to a radio or to an X-ray source for inclusion in this 
sample. However, if an object's optical SDSS position 
matches within 2" to a FIRST radio source, within 10" 
to an NVSS radio source, and/or within 60" to a RASS 
X-ray source, then we keep track of its radio and X-ray 
properties. We also note that radio and/or X-ray infor- 
mation can be used for some ser endipity SDSS spectro - 
scopic target algorithms (e.g., see I Anderson et anT 2003). 
and we consult radio/X-ray information to avoid reject- 
ing a small fraction (~2%) of BL Lac candidates with 
optical colors similar to featureless stars (see fc|4.1[) or 
high-redshift elliptical galaxies (see N4.3|) . Thus, multi- 
wavelength properties play a minor role in our selection 
technique, but the flux sensitivity of this sample is dom- 
inated by the SDSS optical flux limit. 

2.2. SDSS Database Queries 

We first query the SDSS spectroscopic database for 
all objects flagged as having unreliable pipeline redshifts 
(i.e., ZSTATUS = ' N0T_MEASURED ' or 'FAILED'). This 
query returned 12,025 spectra, and we refer to it as the 
'badZ' query. Next we search the database for spectra 
with pipeline redshifts not nagged as ' N0T_MEASURED ' or 
'FAILED' . In this query, we retain spectra with pipeline 



4 



Plotkin et al. 



measured C < 0.41 and pipeline measured REW < 10 A 
for the following spectral lines: Lya, C IV, C III], Mg II, 
[O II], H/3, [O III], and Ha. This yielded 149,321 spectra, 
and we refer to it as the 'goodZ' query. Not all pipeline 
rcdshifts returned by this query are correct. However, 
objects with incorrect redshifts are almost always es- 
pecially featureless and are retained. Although SDSS 
pipeline equivalent width measures are typically reliable, 
they can sometimes fail for especially weak features; this 
is the reason we initially require REWs twice as large as 
our adopted BL Lac classification criteria. We only ex- 
tend the H/K break criterion from C < 0.40 to C < 0.41 
during the database queries because the pipeline H/K 
break measures are fairly reliable (even for small H/K 
breaks): independent measurements of the H/K breaks 
of ~200 objects are nearly identical to the pipeline mea- 
sures (<7rms ~ 0.008). 

As expected, these database queries recover the major- 
ity (>90%) but not all previously known SDSS BL Lac 
candidates. For example, 224 of C05's 240 probable 
BL Lac candidates are returned by our database queries; 
simil arly, 465 of the 501 ra dio-selected BL Lac candidates 
from lPlotkin et all (|2008f ) are recovered. Further details 
are provided in §5.11 

2.3. S/N Cut 

We remove especially noisy spectra from our sample 
by requiring each spectrum to have a signal-to-noise ra- 
tio S/N > 100 in at least one of 3 spectral regimes cen- 
tered on 4750, 6250, and 7750 A with width AA = 500 A. 
Each wavelength region typically contains approximately 
450, 350, and 300 pixels respectively. These spectral re- 
gions are typically of high quality and fall within the g, r, 
and i filters respectively, and this S/N constraint roughly 
corresponds to fiber magnitudes g < 20.5, r < 20.3, or 
i < 19.6 for a typical spectroscopic plate. These cuts re- 
duce the number of spectra returned by the 'badZ' and 
'goodZ' queries to 2197 and 140,493 spectra respectively. 
More sources are eliminated from the 'badZ' query be- 
cause spectra for which the SDSS pipeline does not find 
reliable redshifts tend to have poorer S/N. 

2.4. Automated Removal of the Majority of Stars and 

Galaxies 

We are now left with 142,690 spectra, of which, given 
the area of the sky surveyed, we expect <0.5% to be 
BL Lac objects. The bulk of the returned objects are 
expected to be stars, especially because about 10% of 
the DR7.1 spectra are taken for SEGUE. 13 Our queries 
also return quasars with weak emission features, galaxies 
that underwent recent star formation, and some passive 
ellipticals that have exceptionally weak H/K breaks. It 
is impractical to uniformly visually inspect this many 
spectra, so we seek additional automated algorithms to 
reduce the number of spectra that will require manual 
inspection. 

Following |C0l (see their §3.1.2), we make use of the 
'SpecBS' SDSS spectral classification routines. 14 The 
'SpecBS' routines fit each spectrum with various stellar 

13 It is possible for BL Lac objects to appear on SEGUE plates, 
which is why we search these plates despite the increased stellar 
contamination. 

14 The 'SpecBS' routines were written by David J. Schlegel. 



templates (to classify stars) and combinations of eigen- 
spectra (to classify galaxies and quasars). A second or 
third-order polynomial is added to the fits to account 
for deviations of individual spectra from the stellar tem- 
plates and eigenspectra. The spectrum's redshift is left as 
a free parameter, and the stellar templatc/cigcnspcctra 
plus polynomial fit that yields the lowest % 2 is deter- 
mined as the best fit. The source is classified as a star, 
as a galaxy, or as a quasar according to the preferred 
stellar template or combination of eigenspectra. 

We use the strength of the polynomial component to 
remove contaminating spectra. Normal stars (i.e., with 
modest absorption features) and normal galaxies should 
be well fit by stellar templates and galaxy eigenspectra 
with no need for a strong polynomial component. Since 
the odd spectral characteristics of BL Lac objects are 
caused by non-thermal radiation, BL Lac spectra that 
are classified as stars or galaxies by the 'specBS' routines 
will require significant polynomial contributions. There 
are combinations of quasar eigenspectra that can yield 
non-thermal spectra with weak emission features with- 
out strong polynomial components, so the strength of 
the requisite polynomial cannot be used to reject spec- 
tra classified as quasars by the 'specBS' routines. How- 
ever, only ~10% of all SDSS spectra are classified as 
quasars, so applying a cut based on polynomial strength 
to even just the objects classified as stars and galaxies 
will greatly reduce the number of spectra to be examined 
manually. 

Following [C0l we calculate the relative strength of 
the polynomial component in two spectral regions. The 
first region is near the blue end of the SDSS spectrum 
(3800-4200 A), and the second is near the red end (8800- 
9200 A) . For each spectral regime we measure the average 
total flux density, {ftot), and we calculate the average flux 
density from the polynomial, (f po iy)- We then define the 
polynomial ratio R = (f po iy) / {ftot)- We reject all ob- 
jects classified as stars or galaxies for which Ruue < 0.6 
and R re d < 0.5 (i.e., we reject objects if the polynomial 
does not account for at least half of the observed flux). 
The higher cutoff in the blue is due to the larger number 
of absorption features in stellar spectra and the pres- 
ence of the H/K break for some galaxies in this spectral 
regime. As a consistency check, we visually inspected a 
random subset of spectra that failed the above polyno- 
mial ratio cut, and those spectra indeed appear to be 
normal galaxies and stars. 

We do not perform this polynomial test on the 2197 
remaining spectra from the 'badZ' query; we prefer to 
examine all of those visually since that might lead to 
serendipitous discoveries of other classes of objects. We 
also retain all objects returned by the 'goodZ' query that 
are classified as quasars (2832 spectra) by the 'SpecBS' 
routines. We apply the polynomial cut to objects classi- 
fied as stars (81,883 spectra) and galaxies (55,778 spec- 
tra). This cut retains 21,049 spectra from the 'goodZ' 
query (2,832 classified as quasars by 'SpecBS', 14,909 
stars, and 3,308 galaxies). The selection process until 
visual inspection is summarized in Figure [T] 

A brief description is given in §4.2 of Adelman-McCarthv et al. 
(2007), and starting with DR6 the 'SpecBS' data products are 
available for download from the SDSS Data Archive Server (DAS, 
http : //www. sdss . org/dr7/access/index .html#DAS ). 
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Fig. 1. — Flow chart s umm a rizin g the selection process until 
visual inspection (Sections 12.21 — f2. 41 1. 

3. VISUAL INSPECTION 

We are left with 23,246 spectra to examine visually. 
The vast majority of these spectra are stars, and most of 
them could be rejected through automated routines that 
search for spectra with common stellar absorption fea- 
tures. However, stars that would be rejected with such 
an algorithm can also be very easily and quickly identi- 
fied by eye. We thus opted to visually inspect all ^23,000 
spectra to expand the opportunity for potential serendip- 
itous discoveries. We also examine color cutouts of the 
SDSS images 15 to reject objects that are obvious super- 
positions, blends, or faint low-surface brightness galaxies 
(which can have weak featured spectra that formally pass 
our BL Lac selection criteria). 

During visual inspection we attempt to identify a red- 
shift for each spectrum, and we then proceed to classify 
objects as BL Lac candidates, stars, galaxies, or quasars. 
Because we do not consult all SDSS spectral quality flags 
during our database queries (i )2.2p . a small percentage of 
spectra (~l-2%) are missing data over portions of the 
SDSS wavelength coverage. We reject the small fraction 
of spectra that are missing substantial portions of data, 
and we reject them outright (in addition to the spec- 
tra we reject based on their S/N). Our visual inspection 
procedure is outlined below and summarized in Figure [2j 

3.1. Spectroscopic Redshifts 

We verify the pipeline redshift by searching for at least 
two spectral features consistent with the pipeline red- 
shift. If the pipeline redshift is correct, or if it is incor- 
rect but we identify a different redshift consistent with at 
least two spectral features, then wc refer to the redshift 
as reliable. Spectra are also assigned reliable redshifts if 
the Lya forest is the only observed spectral feature. If 
the pipeline redshift (or our best estimate) appears cor- 
rect but the features are especially weak, or the spectrum 
is too noisy to claim the redshift with high-confidence, 
then we classify the redshift as tentative. A spectrum can 
also be assigned a tentative redshift if it shows only a sin- 
gle weak emission feature that we assume to be Mg II. 
A spectrum that shows an absorption doublet consistent 
with Mg II A2796,2804 (that might be intrinsic to the 
source, or caused by intervening material) is said to have 

15 http: //cas . sdss . org/astrodr7/en/tools/chart/list . asp 
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Fig. 2. — Flow chart summarizing the visual inspection process 

a lower limit redshift. Note, we require detection of both 
Mg II components, so we are not sensitive to interven- 
ing absorption systems with broad (i.e., blended) Mg II 
absorption. Finally, a spectrum that is essentially fea- 
tureless or too noisy to determine a redshift with any 
confidence is said to have an unknown redshift. Objects 
with unknown redshift are constrained to have z < 2.2, or 
else they would show the Lya forest in their SDSS spec- 
tra. In the following, we treat all tentative and lower 
limit redshifts as exact. 

3.2. Spectral Line Measurements 

We measure C (the strength of the Ca II H/K break, 
see Equation [TJ via an automated process for all spectra 
with reliable, tentative, and lower limit redshifts; spec- 
tra with H/K breaks larger than 40% are rejected. The 
REW of the strongest emission line in each spectrum is 
measured manually with IRAF, and spectra showing any 
emission line with REW > 5 A are rejected. We fit se- 
lected blended lines with multiple Gaussians. However, 
we reject spectra that have deblended Ha+[N II] compo- 
nents with REW < 5 A if they show blended Ha+[N II] 
with REW > 5 A and broad Ha. This is under the 
assumption that deblending underestimates the Ha line 
flux and overestimates the [N II] line flux. 

3.3. Contaminants Removed During Visual Inspection 
3.3.1. Stars 

We reject spectra that show obvious stellar absorption 
features at zero redshift (e.g., the Balmer series or He for 
O, B and A stars; Ca II, Mg I, and Na I for F, G and K 
stars; Ca II, Na I, and TiO molecular bands for M stars, 
etc.). Approximately 70-75% of spectra are rejected as 



6 



Plotkin et al. 




4000 6000 8000 4000 6000 8000 

Wavelength (A) 



Fig. 3. — Sample spectra rejected during visual inspection. Se- 
lected spectral features are labeled; absorption features are marked 
with dashed lines and emission features are marked with dotted 
lines. Top left: SDSS J080131. 15+532900.8, a DZ white dwarf 
with prominent Ca II absorption but an otherwise relatively fea- 
tureless spectrum. Top right: SDSS J122801. 59+330035. 5, a mod- 
erately cool DC white dwarf candidate with a relatively featureless 
spectrum and a blackbody continuum that peaks in the SDSS spec- 
tral coverage. Bottom left: SDSS J113338.64+041146.2, an E+A 
galaxy rejected because it shows H<5 with REW < —5 A in absorp- 
tion. Bottom right: SDSS J144408.82+585332.1, a quasar with 
Mg II emission cut by absorption. 

obvious stars during visual inspection. The large stellar 
contamination rate illustrates one of the disadvantages 
of not using radio and/or X-ray emission as selection 
criteria. Less than 10% of spectra that r equired visual 
inspec tion in the radio selected sample of iPlotkin et alJ 
(2008) were identified as stars; these star s are likely su- 
perpo sitions with radio sources (e.g., see iKimball et al.1 
2009). Of course, inclusion of such multiwavelength se- 
lection criteria likely biases samples toward radio and / or 
X-ray bright objects. In any case, the removal of even 
~10 obvious stars is not prohibitively time consuming. 

Main sequence stars of all spectral types were iden- 
tified as contaminants. We also identified about a 
dozen cataclysmic variables and a similar number of M- 
dwarf/ white dwarf binaries, all of which were previously 
known. Thousands of white dwarfs were additionally re- 
covered, including some members of unusual white dwarf 
subclasses: these include 8 DZ white dwarfs with promi- 
nent Ca II absorption (but otherwise relatively feature- 
less spectra) that are not identified as white dwarfs in 
SIMBAD; 16 we also recovered relatively featureless stars 
with blackbody continua that peak in the SDSS spec- 
tral coverage. These are likely cool DC white dwarfs; 
we find 9 that are not cataloged as such in SIMBAD. 
The coordinates of these 17 serendipitous discoveries are 
listed in Table [TJ and sample spectra of two white dwarf 
candidates are shown in the top row of Figure [3] 

We fit blackbodies to the spectra of these 17 white 
dwarfs (excluding Ca II regions for objects with strong 
absorption), and we obtain good fits to each star. The 
best fit temperature estimates are listed in Table [TJ Note 
that none of the cool white dwarfs are at low enough 
temperature to be confidently considered "ultra-cool" 
(T < 4000 K), which is also evidenced by the fact that 

16 The SIMBAD database (http://simbad.u-strasbg.fr/simbad/ ) 

is operated at CDS, Strasbourg, France. 



their continua are well-fit by simple blackbodies. Ultra- 
cool white dwarfs have significant sources of opacity in 
their atmospheres from collisions between H2 molecules 
(or between H2 and He), and simple blackbodies there- 
fore do not fit their spectra well w ithout detail ed model- 
ing of their atmospheres (e.g., see lFarihill2005h . 

Visual inspection does not remove contamination from 
especially featureless stars like DC white dwarfs too hot 
for their continua to peak in the SDSS spectral coverage. 
Their removal is discussed in §4.11 

3.3.2. Galaxies 

A spectrum is rejected as a galaxy for any of the fol- 
lowing reasons (assuming it has a reliable, tentative, or 
lower limit redshift): 

1. its Ca II H/K break is larger than 40%; 

2. any emission feature has REW > 5 A; 

3. it shows H<5 in absorption with REW < —5 A. 
This criterion is designed to remove galaxies that 
recently underwent a period of star formation but 
no longer sho w strong emission features (e.g., E+A 
galaxies; see iGoto et al.1 12003). A sample E+A 
galaxy spectrum rejected based on H<5 is shown in 
the bottom left panel of Figure [3J 

4. the SDSS cutout image shows an obvious blue, low 
surface brightness galaxy. These galaxies, espe- 
cially because they tend to have poor S/N, can 
possess featureless spectra within our quantitative 
limits. A sample spectrum and SDSS image of a 
rejected low surface brightness galaxy are shown in 
the top left panel of Figure |U 

With the above recipe, we reject 5-10% of the visually 
inspected spectra as galaxies. Some passive elliptical 
galaxies with exceptionally weak H/K breaks might still 
contaminate our sample. However, any such population 
should be extremely small. 

The above recipe does not reject a small population 
of ~80 faint and extremely red galaxies (g — r > 1.4) 
at z > 0.2. These galaxies were targeted for SDSS spec- 
tro scopy by the Luminous Red Galaxy algorithms (LRG, 
see lEisenstein et all 120011 ); their removal is discussed in 

S3 

3.3.3. Quasars 

Quasars are rejected if any standard emission feature 
(e.g., Lya, C IV, C III], Mg II, H/3, or Ha) has REW > 
5 A. This may, however, fail to reject quasars that are 
weak in such standard emission lines but show unusually 
strong blends of iron emission features. We also reject 
objects with obvious broad absorption features, or for 
which we measure REW < 5 A but the emission feature 
is clearly cut by absorption (see the bottom right panel of 
Figure [3]). Approximately 10% of the visually inspected 
spectra are rejected as quasars. 

3.3.4. A Strange Class of BALQSOs 

We serendipitously recovered/discovered 9 higher- 
rcdshift (1.1 < z < 1.9) objects that show extreme drop- 
offs in their continua blueward o f rest-fram e 2800 A, simi- 
lar to two objects discovered bv lHall et al.1 (|2002f ). These 
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TABLE 1 

Unusual Subclasses of White Dwarfs 



SDSS Name 




RA 


Dec 


r 




Temperature 


Comment 




(J2000) 




(J2000) 


(J2000) 


(ma 


?) 


(K) 






080131.15+532900 


8 


120.37982 


53.48358 


17 


19 


9719 


Strong Ca 


11 


Uo2927.o5+U75911 


1 


127.OO0U5 


■7 n r c 

7. 98052 


18 


69 


11232 


Strong Ca 


11 
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1 


1 07 one 07 


20.2yo24 


20 
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9369 
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UoOOOO.OO + 2o22o2 
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1 on 1 a a m 
lzy. f 44U1 


no 001 00 
ZO.OOlOO 
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86 


6625 
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11 


UyUz4U.4 7 +I0000D 
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1 or c^O^ A 

1do.oooo4 


1 k c:nm 1 
15.59911 


17 


65 


4467 


Cool DO 




Uy24oU.oU-|-olZUo2 
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1/11 1 oqq c; 
141. IzooO 


Q1 Q/IOQI 

ol .o4zol 
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40Do 


OOOl DO 




094911.91+182231 


4 


147.29966 


18.37541 


19 


98 


10157 


Strong Ca 


11 


111109.10+384856 


8 


167.78795 


38.81578 


18 


21 


5298 


Cool DC 




111316.48+285905 


6 


168.31869 


28.98492 


17 


72 


4726 


Cool DC 




120435.37+231607 


2 


181.14739 


23.26868 


18 


58 


4960 


Cool DC 




121051.65+313659 


9 


182.71523 


31.61664 


18 


93 


9841 


Strong Ca 


11 


122801.59+330035 


5 


187.00664 


33.00988 


17 


99 


4483 


Cool DC 




124805.38+141140 


3 


192.02245 


14.19455 


19 


59 


13345 


Strong Ca 


11 


125307.45+220359 


1 


193.28107 


22.06644 


18 


10 


5226 


Cool DC 




133319.26+245046 


6 


203.33029 


24.84628 


18 


57 


4434 


Cool DC 




134156.62+133845 


3 


205.48595 


13.64592 


19 


24 


11998 


Strong Ca 


11 


162356.37+133614 


1 


245.98490 


13.60392 


18 


66 


4626 


Cool DC 
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Wavelength (A) 

Fig. 4. — More sample rejected spectra, with selected spectral ab- 
sorption features labeled. Corresponding SDSS color image cutouts 
are shown for the top 2 spectra, as these images were consulted for 
their removal during visual inspection. The SDSS images are 50" 
along each side. Top left: SDSS J111701. 05+130558.7, rejected 
during visual inspection as a nearby low-surface brightness galaxy. 
Top right: SDSS J152036.48+133407.3, rejected during visual in- 
spection as a superposition of a (red) galaxy and a (blue) star. 
Both components fall within the 3" SDSS optical fiber. Selected 
stellar absorption features are labeled and marked with dashed 
lines. Selected absorption features from the galaxy are labeled 
and marked with dotted lines. The galaxy features are all consis- 
tent with z = 0.16. Bottom left: SDSS J212745. 21-074359.6, re- 
jected post-visual inspection as an LRG (see JO), Bottom right: 
SDSS J082637. 85+503357. 3, rejected post-visual inspecti on a s a 
DC white dwarf candidate based on its optical colors (see £|4.1| |. 

objects are examples of reddened and/or overlapping- 
trough broad absorption line quasars (BALQSOs). A 
similar object (VPMS J1342+2840; z ~ 1.3) that does 
not appear in the SDSS spe c trosco pic database was dis- 
covered by iMeusinger et all (|2005f ). and two similar low 
rcdshift analogs (FBQS 1503+2330; z=0.49 2 and FBQS 
1055+ 3124; z=0.404) are presented in iWhite et all 
(2000). All previously known objects in this class are ra- 
dio emitters, but they lack X-ray detections at the RASS 
sensitivity level. None of the 9 objects presented here 
have X-ray detections in RASS, and all but two have ra- 
dio detections in FIRST. The two objects lacking FIRST 
detections are not detected in NVSS either. 



Basic properties of the 9 recovered objects are given in 
Table [2 and their spectra are shown in Figure [5l The 
sources are presented in order of decreasing strength of 
Mg II absorption. Included in these 9 objects (for com- 
pleteness) are the two prototypes ( SPSS J0105— 0033 and 
SDSS J2204+0031) discovered bv lHall et all (|2002T ). as 
well as two similar objects (SDSS J 1309+1125 and SPS S 
J1450+4615) initially presented in lPlotkin et all (|2008f) . 
We also note in Table [2] whether these objects are cat- 
aloged as quasars in SPSS quasar catalogs, or if they 
appear in SPSS BALQSO catalogs. (An object's previ- 
ous appearance in an SPSS catalog, however, may not 
indicate that it has been pre viously associat ed with the 
strange objects discovered bv lHall eta l. 2002). Each red- 
shift is derived from a Mg II absorption doublet, except 
for SPSS J1618+1533. We assign this object a redshift 
z = 1.359 from [O II] emission and its continuum shape; 
we interpret Fe II and Mg II absorption in its spectrum 
to be from an intervening system at z = 1.156. 

The increased sample size of these rare and unusual 
objects, now covering a broader redshift range, may allow 
for tighter constraints on their true nature. We do not 
include any of these 9 objects as BL Lac candidates. 

3.3.5. Weak Line Quasars 

The SPSS has discovered a population of ~80 high- 
redshift quasars with weak emission lines (quantified as 
Ly q+N V with BMW < 10-15 A) at 2.7 < z < 5.9 (e.g. 
seelFan et al.lll999l:lAn"derson et a ll200lH She mmer et all 
l200fll200aiDiamond-Stanic et al.ll2009f) . Their redshifts 
can typically only be derived from the onset of the Lya 
forest. While their optical spectra are reminiscent of 
BL Lac objects, WLQs tend to be weaker radio and X- 
ray emitters, and they do not show strong polarization. 
There does not appear to be a large parent population of 
radio-loud BL Lac objects at such high-redshifts, and it is 
unlikely that the SPSS is biased against recovering high- 
rcdshift radio and X-ray bright o bjects with featureless 
optic al sp ectra (e.g., see §4.1 of IShemmer et "all l2?M 
and i j6.2l of this paper). Thus, WLQs' lack of strong 
spectral features might not be caused by relativistically 
boosted jet emission. However, because WLQs formally 
pass our BL Lac criteria, we include all high-redshift ob- 
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TABLE 2 

Unusual Broad Absorption Line Quasars 



SDSS Name 


RA 


Dec 


Redshift 


r 


Radio Flux a 


Refer 


(J2000) 


(J2000) 


(J2000) 




(mag) 


(mJy) 




155626.90+103457.3 


239.11212 


10.58260 


1.889 


19.45 






134951.94+382334.1 


207.46642 


38.39281 


1.094 


19.13 


1.73 


3; 5 


160827.08+075811.5 


242.11286 


7.96988 


1.183 


16.94 


13.79 


3 


145045.56+461504.2 


222.68984 


46.25118 


1.877 


19.54 


1.51 


3;4 


220445.26+003141.9 


331.18859 


0.52832 


1.335 


17.35 


3.03 


1;2;3 


130941.36+112540.1 


197.42234 


11.42783 


1.363 


18.57 


0.83 


3;4;5 


010540.75-003314.0 


16.41980 


-0.55389 


1.179 


17.97 


4.84 


1;2;3 


215950.30+124718.4 


329.95960 


12.78846 


1.514 


21.03 




2;3;5 


161836.09+153313.6 


244.65041 


15.55378 


1.359 b 


18.91 


2.10 





Reference s. — (1) Identified as an unusual BALQSO in IHall et all J2002|): (2) cataloged as a 
BALQSO in ITrump et al.l d2006D : (3) catalo ged as an SDSS qua sar in ISchneider et al] d2007D : (4) 
identified as an un usual BALQSO quasar in Plotkin ct al. (2008); (5) cataloged as a BALQSO in 
IGibson et"aH 1 120091) . 

a Integrated Flux Density at 1.4 GHz from FIRST. 

Also shows intervening Mg II and Fe II absorption at z — 1.156. 



E 
o 

cn 



CO 

c 
CD 

Q 



6 


\ ^ 




T — 1 — 1 — ■ — 1 — ■ — 1 — ' 


1 1 ' 


■ — i — ■ — 

50 


^ — i — i — ■ — i — i — ■ — ' 


1 1 1 ' 1 "j 


4 




s : 






40 


■ J 


X | 


2 










VM 


wu 

































i u 




20 






-2 






2 

z = 1.89 




z = 


= 1.09'° 




z = 


1.18 1 


6 
























r\ 








4 




30 




















\J 




10 






2 






20 








| fj 





























10 
















z = 1.88 




z = 


= 1.34 




z = 


1,36 : 






10 


;i ■ ■ ■ 1 ■ ■ ■ 


H— 










20 












8 


























15 


\ J 






1 


6 


: J 


\ 




10 










4 












2 












5 










2 








J 1 1 1 1 1 1 L 


z = 1.18 o 

_l 1 1 I 1 1 : 


J D 1 1 1 1 1 t 


z = 
i i 




.51 

. i . 


j D 1 1 1 ■ ■ X 


z = 
1 1 1_ 


1.36 : 



4000 1000 2000 3000 .4000 1000 2000 3000 



Wavelength (A) 

Fig. 5. — SDSS spectra of reddened and/or overlapping-trough BALQSOs, ordered by decreasing strength of Mg II absorption, 
spectra are shown in their rest-frame, and the location of Mg II is marked with a vertical dotted line in each panel. 



The 



jects with REW < 5 A without prejudice (but we do 
identify a subset of objects with weak radio-emission in 

Part of the motivation for assembling an optically se- 
lected BL Lac sample is to find interesting objects like 
WLQs that would not be recovered by radio and/or 
X-ray selection. In Table [3] we list 20 serendipitously 
discovered WLQs (that have reliable redshifts and are 
ne w SDSS WLQ identification s, i. e., they do not appear 
in iShemmer et all [2006, 20Ql or iDiamond-Stanic et al.1 
12009ft with REW < 10 A for Lya+N V. For 6 of these 
objects, all emission features show REW < 5 A, and 
we therefore include those 6 objects as BL Lac candi- 
dates as well. The remaining 14 objects are not included 
as BL Lac candidates. All 20 objects in Table [3] have 
reliable z > 2.2, which is where the Lya forest enters 



the SDSS spectral coverage. We also note that the high- 
redshift nature of WLQs discovered by the SDSS is a 
selection effect because the Lya forest is used to de- 
termin e redshifts. For example, IDiamond-Stanic et al.l 
(|2009f) only searched for WLQs at z > 3. There are 
some lower- redshift radio- weak/quiet objects with weak- 
featured s pectra in the literatu re: e.g., PG 1407+265 
(z = 0.94. iMcDowell et aTl ll995); PHL 1811 (z = 0.192, 
Leighlv et al.ll2007n . 2QZJ215454.3-305654 \z = 0.494, 



Londish et al.n20?a: S PSS J094533.99+100950.1 {z 
1.66, ilr yniewicz et al.l l2~009[ ) . Our BL Lac selection al- 



g orithm is also sensi tive to such objects (see _ 

iSmith et all (|2007t ) found that optically selected SDSS 
BL Lac candidates with z > 1 do not show strong polar- 
ization (< 3%), suggesting that some z > 1 BL Lac can- 
didates might not have continua dominated by beamed 
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TABLE 3 

New z > 2.2 Weak Line Quasar Identifications 



SDSS Name 


RA 


Dec 


Redshift 


z 






Radio Flux° 


a ro d 


(J2000) 


(J2000) 


(J2000) 




(mag) 


(mag) 




(mJy) 




080906.88+172955.2° 


122.27867 


17.49867 


2.951 


18.17 


0.07 


0.47 


<0.92 


<0.09 


082638.59+515233.2° 


126.66080 


51.87591 


2.850 


16.93 


0.05 


0.67 


1.74 


0.04 


083304.74+415331.3° 


128.26976 


41.89204 


2.329 


18.92 


0.05 


0.78 


<0.95 


<0.15 


083330.56+233909.2 


128.37736 


on f* f o f* 

23.65256 


2.417 


18.90 


0.05 


1.00 


2.05 


0.20 


Uo4z49.Uz+zo5zU4.7 


ion 7n a on 

130.70429 


oo o/"onn 

23. 808OO 


3.316 


18.98 


0.04 


1.46 


1.55 


0.12 


084424.24+124546. 5 f 


131.10102 


12.76294 


2.466 


17.60 


0.07 


1.03 


4.17 


0.14 


nr\n' h mo no i .4 "i ot a o oo 

090703. 92+410748. 3 L 


13d.7do34 


a "i "i o nnn 

41.13009 


2.672 


19.33 


0.02 


0.08 


<1.39 


<0.25 


Oyzolz. / o+l / 4452. 5 


140.o0ol4 


17 7< on 1 


2.260 


18.21 


0.05 


1.11 


<0.95 


<0.08 


nno^oT tr o i o c o o o o j^i 

093437. 53+2o2232. 6 


143.65639 


O/" 1 OTC T A 

2b. 37574 


3.370 


19.92 


0.03 


0.76 


<0.90 


<0.19 


101849.78+271914.9 


154.70745 


27.32081 


2.603 


18.95 


0.04 


1.09 


<1.04 


<0.13 


102609.92+253651.2 


156.54136 


25.61425 


2.318 


18.75 


0.03 


0.54 


3.93 


0.27 


111642.81+420324.9° 


169.17840 


42.05694 


2.526 


18.41 


0.03 


0.80 


1.13 


0.12 


113747.64+391941. 5°- f 


174.44852 


39.32820 


2.395 


18.61 


0.03 


1.12 


<0.90 


<0.09 


115959.71+410152.9° 


179.99880 


41.03137 


2.788 


17.18 


0.02 


0.74 


<0.97 


<0.01 


124745.39+325147.0 f 


191.93914 


32.86308 


2.249 


18.43 


0.02 


0.83 


<0.76 


<0.08 


132703.26+341321.7° 


201.76360 


34.22271 


2.558 


18.62 


0.02 


0.66 


<0.90 


<0.13 


141657.93+123431.6° 


214.24138 


12.57546 


2.609 


18.32 


0.03 


0.53 


<1.00 


<0.12 


144803.36+240704.2 f 


222.01404 


24.11785 


3.544 


18.84 


0.06 


1.03 


<0.94 


<0.08 


213742.25-003912.7° 


324.42605 


-0.65355 


2.257 


19.02 


0.07 


0.93 


<0.91 


<0.15 


233939.48-103539. 3°> f 


354.91454 


-10.59427 


2.757 


18.52 


0.04 


0.41 


<1.02 


<0.14 



a Extinction in the z filter from lSchlegel et al.| J1998D ■ 

Optical spectral index (f u ~ measured via fitting a power law to each SDSS spectrum, excluding Lya 

and the Lya forest from each fit (sec i|5.2t . 
c Integrated Flux Density at 1.4 GHz from FIRST. 

Radio to optical broad-band spectral index referenced at rcst-framcs 5 GHz and 5000 A. We use the SDSS z-band 
magnitude for this calculation because its wavelength coverage is closest to rest-frame 5000 A among the SDSS 
filters. We adopt a ro — 0.2 as the approximate boundary between radio-loud and radio-quiet quasars. Sec SI5.5I 
for details on the a ro calculation. 

° Included in the SDSS DR5 Quasar Catalog jSchneider et al.ll2007l) . 

^ Included as a weak-featured radio-quiet object in Tables [6] and [8] (sec E|5.6t . 



synchrotron radiation. Perhaps z > 1 BL Lac candidates 
are actually a mix of BL Lac objects and lower-redshift 
WLQs. The resolution of this issue is unclear with extant 
data; additional observations (such as variability and 
polarization monitoring, as well as deeper radio/X-ray 
imaging for objects lacking radio/X-ray detections) are 
needed. Thus, we include all z > 1 objects that formally 
pass our BL Lac criteria as BL Lac candidates, but we 
recognize that some of these objects might rather consti- 
tute an equally fascinating population of lower-redshift 
an alogs to WLQs. This issue is discussed in more detail 
in iEl 

3.3.6. Duplicate Spectra and Superpositions 

A small fraction of SDSS objects (<1%) have repeat 
spectroscopic observations. Multiply-observed sources 
only need to pass our selection criteria at a single epoch 
for inclusion in our sample. If multiple spectra of the 
same object are classified as BL Lacs, then we only in- 
clude one (randomly selected) spectrum in our final sam- 
ple. 

The combination of flux contributed by two distinct 
superposed sources can sometimes act to make the ob- 
served spectrum appear featureless. Around 1-2% of the 
manually inspected spectra are random superpositions. 
Superpositions are identified by sets of spectral features 
that are consistent with two different redshifts; some su- 
perpositions are also identified as two objects that fall 
within the 3" SDSS optical fibers, as judged by exam- 
ining their SDSS images. The vast majority of super- 
positions are foreground stars over galaxies, but we do 



see a small number of galaxy/galaxy, galaxy/quasar, and 
star/quasar pairs. An example of a spectrum of a galaxy 
superposed with a foreground star and a thumbnail of its 
SDSS image is shown in the top right panel of Figure [U 
We constructed our database queries to initially retain 
such superpositions, because otherwise we might reject 
real BL Lac objects that appear as blue point sources at 
the centers of elliptical galaxies. 

4. REMOVAL OF CONTAMINANTS THAT SURVIVE 
VISUAL INSPECTION 

Visual inspection reduced our list of nearly featureless 
spectra to 1310 objects. However, there are still surviv- 
ing contaminants, including DC white dwarfs and high- 
redshift LRGs. Their removal is discussed below and 
summarized in Figure 

4.1. Removing Featureless White Dwarfs 

We use proper moti ons from the SD SS+USNO-B 
proper moti on c atalog (Munn ct al. 2004j) to verify the 
result from IC05I that high proper motion objects with 
nearly featureless spectra tend to populate bluer regions 
of optical color space than BL Lac objects. These con- 
taminants are most likely DC white dwarfs. Figure [7] 
shows optical color-color plots for 1222 surviving objects 
with nearly featureless spectra also in the SDSS+USNO- 
B catalog (|Munn et al.l 120041 ). We require each object 
to match to only a single USNO-B source with a 1" 
error circle centered on the SDSS position. Blue plus 
signs denote objects with measured proper motion fi > 
30 milli-arcsec yr _1 , and red triangles have [i < 30 milli- 
arcsec yr _1 . We choose ri = 30 milli-arcsec yr _1 as the 
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Fig. 6. — Flow chart summarizing the removal of stars and galax- 
ies post- visual inspection (@. 



division between large and small proper motions because 
on ly 0.8% of spect r oscop ically confirmed SDSS quasars 
in iSchneider et al.l (|2007[ ) with exactly one match to a 
USNO-B source have /i > 30 milli-arcsec yr -1 . 

Figure [7] indeed confirms that objects with large proper 
m otion tend to have bluer colors. We adopt the color cut 
of IC05L g — r < 0.35 and r — i < 0.13, to remove nearly 
featureless stars. We refer to this blue color space as 
the gri box, and this cut removes 500 objects. Not all 
stars will have large measured proper motions (and ~7% 
of our nearly featureless spectra do not have proper mo- 
tion measures), so we prefer to reject objects based on 
color rather than measured proper motion. It is possible 
for some BL Lac objects to populate the gri box, so we 
retain spectra with blue colors that have reliable, ten- 
tative, or lower limit redshifts and/or that match to a 
radio and/or X-ray source. We retain 61 such objects as 
BL Lac candidates. 

As stated above, objects inside the gri box with un- 
known redshifts are rejected as likely DC white dwarfs. 
However, there are five objects inside the gri box with 
unknown redshifts that also match to a radio and/or X- 
ray source, and we retain this small number of objects 
as BL Lac candidates because of their multiwavelength 
properties. Approximately 84% of our final BL Lac can- 
didates (see fjSj) match to a FIRST/NVSS radio source or 
to a RASS X-ray source. If a similar percentage of very 
blue BL Lac objects have radio or X-ray emission above 
the relevant multiwavelength survey flux limits, then we 
expect our color cut aimed at removing DC white dwarfs 
to only reject ~1 BL Lac object from our final sample. 

In summary, we remove ~500 nearly featureless spectra 
(probably DC white dwarfs) that: 

1. fall within the gri box; 

2. have unknown redshifts; 

3. and do not match to a radio and/or X-ray source. 

4.2. Removing Stellar Contaminants Outside the gri 

Box 

A very small number of objects outside the gri box 
have measured \i > 30 milli-arcsec yr -1 . We reject those 
objects if they: 



1. have fj, > 30 milli-arcsec yr -1 and match to exactly 
one source in the SDSS+USNO-B proper motion 
catalog; 

2. have an unknown redshift; 

3. and do not match to a radio and/or X-ray source. 

The above criteria remove only 4 objects that would 
not have been removed via other cuts. This is our only 
cut based solely on proper motion. Since less than 1% 
of spectroscopically confirmed SDSS quasars show mea- 
sured proper motions larger than 30 milli-arcsec yr , 
with this cut we risk losing < 0.04 BL Lac objects. 

4.3. Removing High- Redshift Luminous Red Galaxies 

There is a small population of higher-redshift objects 
(z > 0.2) targeted for spectroscopy as LRGs that for- 
mally pass our BL Lac criteria after visual inspection. 
These objects occupy a distinct region of optical color 
space, and we separate them phenomenologically as hav- 
ing g — r > 1.4. This division is drawn as a dashed 
line on the relevant panels in Figure [3 where their sep- 
aration from the rest of our nearly featureless spectra 
can be seen. These objects formally pass our S/N cut, 
but they tend to do so only at the redder wavelength 
regions that we consider in W2.3\ This is because the 
Ca II H/K break crosses from the g to the r filter near 
z ~ 0.4. These spectra have acceptable S/N redward 
of rest-frame 4000 A, but they have relatively low flux 
blueward of rest-frame 4000 A; their H/K break measure- 
ments are therefore relatively noisy. For example, the av- 
erage measured uncertainty for the H/K break for these 
nearly featureless objects with g — r > 1.4 is ac ~ 0.04, 
compared to ac ~ 0.01 for all nearly featureless spec- 
tra with g — r < 1.4. We believe the majority of these 
objects are normal LRGs (i.e., with C ~ 0.5), but be- 
cause of their noisier H/K break measures, we find an 
unexpectedly large number of them to have C < 0.4. 

We take precautions to not remove red BL Lac objects 
with g — r > 1.4 from our sample if they exist. We thus 
only remove objects as high-redshift LRGs if they show 
all of the following: 

1. g - r > 1.4; 

2. their SDSS image morphology is extended (i.e., the 
SDSS morphology parameter TYPE = 3); 

3. and they do not match to a radio and/or X-ray 
source. 

We exclude 83 spectra that fulfill the above requirements, 
leaving 17 objects with g — r > 1.4 that remain as BL Lac 
candidates. We visually examined the SDSS spectra and 
images of these 83 rejected spectra, and we confirm they 
indeed look like low S/N, extended LRGs. 

5. THE OPTICALLY SELECTED SAMPLE 

The final sample includes 723 objects. We subdi- 
vide these into 637 radio-loud BL Lac candidates and 
86 radio-quiet objects with nearly featureless spectra 
(defined by having a ro < 0.2, see N5.5P . The radio- 
quiet objects may be a mix of BL Lac objects on the 
radio-faint tail of the BL Lac distribution and a dis- 
tinct alternate class (see discussion below) of exotic AGN 
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Fig. 7. — Optical colors for 1222 nearly featureless spectra with exactly one match in the proper motion catalog and surviving visual 
inspection. Objects with fi > 30 milli-arcsec yr — 1 are denoted as blue plus signs, objects with fx < 30 milli-arcsec yr _1 are shown as red 
triangles. Dashed lines on relevant panels indicate color cuts applied to remove DC white dwarfs (WD), (g — r, r — i) < (0.35,0.13), and 
high-redshift Luminous Red Galaxies (LRG), g — r > 1.4. (a) g — r vs. u — g; (b) r — i vs. g — r; (c) i — z vs. r — i\ (d) i vs. g — r. All panels 
use psf magnitudes not corrected for Galactic extinction. Not shown are ~90 nearly featureless spectra lacking proper motion information 
(or with more than one match in the proper motion catalog.) 



with extremely featureless spectra. We further subdivide 
each subset into high-confidence ('H') and low-confidence 
('L'). We classify an object as low-confidence if it is un- 
clear where to define the continuum near an emission 
feature, and that emission feature can have measured 
REW larger or smaller than 5 A depending on the con- 
tinuum used. We recover 613 and 24 radio-loud high- 
and low-confidence BL Lac candidates respectively; and 
we find 56 and 30 high- and low-confidence radio-quiet 
objects with nearly featureless spectra. Sample spec- 
tra for 4 (radio-loud) BL Lac candidates are shown in 
Figure [5J and basic sample demographics are summa- 
rized in Table |U We compare the optically selected sam- 
ple to previously pub lished SDSS BL Lac ca talogs (i.e., 
Anderson etal] 120031 : lC05t lAnderson et al.l 120071 : and 
Plotkin et al. 1 12008). and we also examine each source's 
classification in NED. 17 From this, we estimate that 
~35% of the optically selected BL Lac candidates are 
new BL Lac identifications, and ~75% of all of our opti- 
cally selected BL Lac candidates were initially discovered 
by the SDSS (either by this or previous studies). 

We loosely refer to the radio-quiet subset as BL Lac 
candidates until they are discussed in more detail in Sec- 
tions 15.51 and 16. 2\ where we compare them to WLQs. 
Distinguishing these objects from WLQs is not trivial, 
especially because 13 of the 86 radio-quiet BL Lac can- 
didates have z spec > 2.2 and pass the criteria to alterna- 

17 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 



TABLE 4 

Optically Selected Sample Demographics 





Radio-Loud 


Radio-Quiet 




H a 


L b 


H a 


L b 


Total Number 


613 


24 


56 


30 


Reliable z Bpec 


216 


15 


22 


16 


Tentative z S pec 


85 


9 


22 


14 


Lower Limit z apec 


42 





7 





Unknown z spec 


270 





5 





Radio Matches 


571 


19 


9 





X-ray Matches 


281 


8 


2 


3 


LBL C 


188 


4 








HBL d 


425 


20 


56 


30 



a High-confidence BL Lac candidate 
Low-confidence BL Lac candidate 

c Low-energy peaked BL Lac objects 
High-energy peaked BL Lac objects 



tively be classified as WLQs. Note, 6 of those 13 radio- 
quiet objects also appear in Table [3] as new WLQ discov- 
eries. The other 14 new WLQ identifications in Table [3] 
have emission lines too strong to pass our stricter BL Lac 
REW < 5 A criterion (and therefore do not appear in 
the optically selected BL Lac sample). 

5.1. Tests on the Completeness of Our Recovery 
Algorithm 

We test our recovery algorithm by correlating our sam- 
ple to the 226 probable BL Lac candidates recovered by 
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Fig. 8. — Sample spectra of 4 BL Lac candidates, with selected 
spectral features labeled. Absorption features are marked with 
dashed lines, and emission features are marked with dotted lines. 
All 4 objects match to FIRST/NVSS radio sources and RASS X- 
ray sources. The spec tral fits from the host galaxy decomposi- 
tion described in £|5.2I are overplotted as red lines. If there is 
significant host galaxy contamination, then the dotted red lines 
illustrate just the power law flux from the AGN. Top left: SDSS 
J092915. 43+501336.1, a high-confidence BL Lac candidate with a 
spectrum too featureless to derive a spectroscopic redshift. Top 
right: SDSS J084712.93+113350.2, a high-confidence BL Lac can- 
didate with reliable z = 0.198 derived from host galaxy spectral fea- 
tures. Bottom left: SDSS J124312. 73+362743.9, a high-confidence 
BL Lac candidate with z > 0.485 derived from intervening Mg II 
absorption. Bottom right: SDSS J164419. 97+454644.3, a low- 
confidence BL Lac candidate with z = 0.225 derived from host 
galaxy spectral features. This object is classified as low- confidence 
because it shows broad Ha that can have measured REW > 5 A 
depending on the continuum adopted during the measurement. 

IC05I that also appear in the DR7.1 database 18 (see 3!l2"|) : 
we recover all but 28 of these objects. Essentially all 
of the missing [C05| objects have spectral features near 
the quantitative BL Lac crite ria l imits. Although our re- 
covery approach is similar to IC051 since we use different 
spectral reductions it is not surprising that some "bor- 
derline" objects do not pass the spectral criteria here. 
Among the 28 missing objects, 2 were rejected by our 
database queries; we also rejected 11 objects because 
they showed poor S/N ( W2.'3]i , they failed the polynomial 
strength test f N2.4p . or they were deemed to have poor 
quality spectra/photometry during visual inspection; 10 
objects were rejected during visual inspection because 
they showed emission features with REW > 5 A in their 
DR7 spectral reductions; and the remaining 5 objects 
were rejected as stars (either during visual inspection or 
after applying the color cuts described in []4]). 

We perform a similar quality assurance test against the 
500 objects in the DR7.1 database that appear in the 
iPlotkin et alj (|2008l ) radio-selected sample. 19 Our opti- 
cal sample does not include 95 of those 500 radio-selected 
objects, of which 35 were rejected by our database 
queries. Again, the majority of missing objects tend to 
have spectral features near our quantitative limits, and 
we measure their features to be too strong in the DR7 
spectral reductions (but their features pass our spectral 
criteria in earlier data release spectral reductions). An 

18 14 of IC05l 's probable BL Lac candidates do not appear in the 
DR7.1 database. 

19 1 object does not appear in the DR7.1 database. 



exception is a small population of 5 spectra that we no- 
ticed to be superpositions of two sources in their DR7 
spectra. Of the other 55 missing spectra, 33 showed poor 
S/N ff2~3]). failed the polynomial strength test ($23}, or 
were deemed to have poor quality spectra during visual 
inspection; 21 objects were rejected because they showed 
emission features with REW > 5 A in their DR7 spec- 
tra; and 1 spectrum was classified as a star during visual 
inspection. 

Given the post-DR6 improvements to the spectropho- 
tometric pipeline, the fraction of previously identified 
SDSS BL Lac candidates (>80%) recovered by this new 
optically selected sample is satisfactory. As expected, the 
missing ^20% are primarily "borderline" BL Lac candi- 
dates. 

5.2. Spectral Decomposition: Estimating AGN Spectral 
Indices, Redshift Limits, AGN Component Fluxes, 
and Host Galaxy Fluxes 

We use each object's SDSS spectrum to measure the 
AGN's optical spectral index, assuming the AGN flux 
follows a power law. If the host galaxy contributes sig- 
nificantly to the total observed flux, then we perform a 
spectral decomposition of the host galaxy from the SDSS 
spectrum. This decomposition additionally allows esti- 
mates of optical fluxes (or limits) for the host galaxy and 
AGN separately. 

If the BL Lac object has an unknown redshift, then (by 
definition) its spectrum does not show significant con- 
tamination from the host galaxy. For the 275 BL Lac 
candidates with unknown redshifts, we fit their spectra 
(via Levenberg-Marquardt least-squares regression) with 
the following power law: 

/ A \ 

f\,total (A) = fx J • (2) 

f\,totai(X) refers to the observed flux density per unit 
wavelength, and f\ is the flux density at reference wave- 
length An = 6165 A (the effective wavelength of the SDSS 
r filter). There are typically ~3800 degrees of freedom 
(dof), and we estimate the goodness of this power law 
(pi) fit with the reduced Xdof, P i- 

If a spectrum has a measured reliable, tentative, or 
lower limit redshift, then the BL Lac host galaxy contri- 
bution to the observed spectrum might be important. In 
addition to fitting Equation [2] (which has 2 free parame- 
ters, the spectral index a\ and the flux density normal- 
ization /a q ), we also perform a power law fit including 
the contaminating flux from the host galaxy, which w e 
assume to be a normal elliptical 20 (see lUrrv et al.|[2000f) . 
Thus, we find the best fit to the following function: 

fxjotai(X) = Ci/ A , eH (A, z) + / Ao (±\ X . (3) 

f\,totai(X) is again the observed flux density, and 
fx,eii(^,z) is the flux density from the elliptical galaxy 
(at the observ ed redshift). We use t he elliptical galaxy 
template from Mannucci et al. (2001) to model the host 
galaxy flux density. This fit has 3 free parameters, the el- 
liptical galaxy flux density normalization (C\ ) , the power 

20 Normal in every sense other than its active nucleus. 
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law flux density normalization (f\ ), and the spectral in- 
dex of the active nucleus (a\). We constrain C\ and 
/a to always be positive, and we estimate the good- 
ness of the host galaxy decomposition with the reduced 
Xdof decomp- We exclude the Lya forest from the fits (for 
spectra with z > 2.2). 

Equation [3] produces much better fits than Equation [2] 
for objects with redshifts derived from host galaxy spec- 
tral features. When the host galaxy becomes negligible 
(as is the case for higher-redshift and/or highly beamed 
BL Lac candidates), Equation [3] is identical to Equa- 
tionO We use the spectral index returned by the fit with 
the smaller Xdo/' m or( ier to minimize the requisite num- 
ber of free parameters. We adopt the value a„ = 2 — ct\ 
returned by whichever fit is preferred. 21 The spectral fits 
are plotted as red lines over the sample BL Lac spectra 
in Figured For the cases where the fit including the host 
galaxy decomposition was preferred, the power law spec- 
trum contributed by just the AGN component is shown 
as a dotted red line. 

The distribution of measured spectral indices is shown 
in Figure Hk,. The average spectral index is (av) = 1.15 
(median ct v = 1.06), with a dispersion of 0.69. The 
spread in the ensemble's a u distribution is much larger 
than the typical measurement uncertainties of er Qiy = 
0.01 — 0.03. It is possible that some of the scatter in 
<7 au can also be attributed to using the same host galaxy 
template for each object. This of course only affects the 
BL Lac candidates that preferred the spectral decompo- 
sition over fitting just a power law fit (around 40% of the 
sample). The dispersion in a v for just the objects that 
preferred the power law fit is similarly large. Variabil- 
ity also contributes to the observed scatter in ct v . There 
are 46 BL Lac candidates with multiple epochs of SDSS 
spectroscopy (about two-thirds have two epochs, most 
of the others have three epochs.) We measured a„ for 
each duplicate spectrum (decomposing the host galaxy 
flux when necessary), and we calculated the difference 
between the largest and smallest measured a v for each 
object. The root mean square of the differences in a v 
is ^0.3, which is significant, but unlikely large enough 
to account for the all of the observed dispersion in the 
entire BL Lac sample. Thus, we conclude that the rela- 
tively large dispersion indicates that BL Lac objects do 
not have universal spectral shapes; the large scatter is 
probably due to a combination of different amounts of 
relativistic beaming among objects, variability, and in- 
nate physical differences (such as black hole mass, accre- 
tion rate, etc). 

5.2.1. Estimating Redshift Limits From Host Galaxy 
Contamination in SDSS Spectra 

BL Lac objects live (perhaps exclusively) in giant el- 
liptical galaxies, and their hosts tend to have similar lu- 



minositics within ±1 ma g (iScarpa et al.ll200rtlUrrv et all 
120001) . ISbarufatti et alj (|2005h note that the dispersion 
in host galaxy absolute magnitudes is narrow enough to 
treat BL Lac host galaxies as standard candles, allowing 
the derivation of "imaging" redshifts if the host galaxy 
can be spatially resolved from the AGN. Using a sub- 
set of 64 objects with Hubble Space Telescope (HST) re- 

21 We report a v rather than ct\ for consistency (and easy com- 
parison) with the bulk of the literature. 
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Fig. 9. — AGN component optical spectral indices, and optical, 
radio, and X-ray flux distributions for BL Lac candidates, (a) The 
optical spectral index «„ for the decomposed AGN flux. The width 
of the histogram indicates that the shape of the AGN's optical 
spectrum is not universal among BL Lac objects, (b) Extinction 
corrected optical magnitudes of the decomposed AGN component 
flux in the SDSS filter closest to 5000 A rest-frame, (c) Loga- 
rithm of integrated radio flux densities at 1.4 GHz. When avail- 
able, flux densities are taken from FIRST, otherwise flux densities 
are taken from NVSS. The hatched histogram shows radio flux 
density upper limits from FIRST for objects lacking radio detec- 
tions in both surveys, (d) Logarithm of broad-band X-ray fluxes 
(corrected for Galactic absorption) from 0.1-2.4 keV. The hatched 
histogram shows X-ray flux upper limits for objects lacking X-ray 
detections in RASS. 

solved host galaxies and known spectroscopic redshifts, 
they show that host galaxy apparent magnitudes are 
tightly constrained with redshift on a Hubble diagram 
if they have absolute magnitudes Mr = —22.9 ± 0.5. 22 
Their standard candle algorithm is calibrated out to 
z ~ 0.7, and it provides "imaging" redshifts accurate 
to a z ~ ±0.05 if a host galaxy is resolved; their tech- 
nique can also be used to estimate redshift lower limits 

in the cases of host gala xy non-detections . 

We use the result from Sbarufatti etaLl (f2005h to esti- 
mate redshift lower limits for objects in our SDSS sam- 
ple that lack spectroscopic redshifts. 23 However, instead 
of attempting to spatially resolve host galaxies (which 
becomes more challenging with increasing redshift, espe- 
cially with ground-based telescopes), we rather take ad- 
vantage of each object's high-quality spectrum. With ex- 
tant SDSS data, a spectral decomposition is likely more 
sensitive than a spatial decomposition: spectral decom- 
position is primarily limited by the spectral resolution 
(and collecting area) of the telescope, while spatial de- 
composition limitations are dominated by the spatial res- 
olution. For example, only ~200 SDSS BL Lac candi- 
dates appear extended in their SDSS images, and es- 



22 |Sbarufatti ct al. (2005) use H = 71 km s" 1 Mpc" 
0.7, and Q m = 0.3, similar to our adopted cosmology. 

23 These objects have redshift upper limits z < 2.2, or else they 
would show the Lya forest. 
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sentially all of these objects already have spectroscopic 
redshifts (and they show host galaxy contamination in 
their spectra.) However, some objects with unresolved 
host galaxies in their SDSS images still show host galaxy 
features in their SDSS spectra. 

Our goal is to estimate an i?-band apparent magnitude 
limit (in the Cousins filter system) for each host galaxy. 
Equipped with these flux estimates, we can then trans- 
late apparent magnitudes into redshift limits by assum- 
ing the hosts are standard candles. In the previous sub- 
section, approximately 200 high-confidence BL Lac can- 
didates with reliable redshifts preferred the host galaxy 
decomposition fit (Equation [3]) over just a power law fit 
(Equation [2]) . Over 90% of these objects have decom- 
posed AGN flux densities at most 5 times larger than 
their host galaxy flux den sities at 6581 A (the effective 
wavelength of the R filter: iFukugita et al.lll996T) . There- 
fore, for the 275 BL Lac candidates lacking spectroscopic 
redshifts, we assume their host galaxy flux densities are 
at least 5 times fainter than the observed flux densities 
at 6581 A. Otherwise we would have seen host galaxy 
features in their spectra during visual inspection and as- 
signed them spectroscopic redshifts. 

We place upper limits to the host galaxy flux densities 
per unit frequency at the effective wavelength of the R 
filter (X e ff,R = 6581 A) as: 



fu,hg(z) < 0-2 X 



X 2 



A( z )f\,total{Kff,R.) 



, (4) 



where f\ t totai(X e ff,R) is the observed flux density per 
unit wavelength at X e ff,B- A(z) is an aperture correc- 
tion, to correct for the amount of host galaxy flux that 
does not fall within each SDSS optical fiber (which are 
3" in diameter.) The fraction of light lost due to the 
finite fiber size can be significant, especially at lower 
redshifts where the host galaxy appears more extended. 
We estimate the size of the galaxy as a function of red- 
shift, assuming a de Vaucouleurs profile with effective 
radius r e = lOkpc (typical of BL Lac host galaxies, see 
lUrrv et aT]|2000f ). A(z) is then approximately the total 
galaxy flux integrated over an infinite aperture divided 
by the amount of light expected to fall within 3" at each 
redshift (accounting for cosmological effects.) Note, since 
the AGN is a point source, no aperture correction needs 
to be applied to its flux estimate. 

We convert f v ,hg{z) to a limit on the i?-band host 
galaxy apparent magnitude, Rh g {z), using the flux den- 
sity zero point of the AB magnitude scale (3631 Jy 
at \ e ffM) ■ The redshift for which Rh g (z) best 
matches the standard candle relation in Equation 2 of 
ISbarufatti et a l. (2005) is then taken as a lower limit to 
the objects' redshift, Zh g - We refer to these redshift lim- 
its throughout the text as "host galaxy" redshift limits. 
The standard candle relation is only calibrated out to 
z ~ 0.7, so we assign lower limits of Zh g = 0.7 for 15 
objects whose host galaxy apparent magnitudes are al- 
ways fainter than expected from the standard candle re- 
la tion. Our approa c h is si milar to the technique outlined 
in ISbarufatti et all (|2006f) . who derive redshift limits by 
constraining AGN to host galaxy flux ratios in feature- 
less BL Lac spectra using the equivalent width detection 
threshold of their high-resolution spectra. 



The redshift distribution for all 723 objects (i.e., in- 
cluding all spectroscopic redshifts and "host galaxy" red- 
shift limits) is shown as open histograms in panels (a) 
and (c) of Figure [TO] The distribution of spectroscopic 
redshifts (assuming tentative and lower limit redshifts 
are exact) is shown as the hatched red histogram in 
panel (a). Lower limit spectroscopic redshifts (49 ob- 
jects) are shown as the solid red histogram in that panel, 
and the inset shows a blowup of z spec < 0.7. The dis- 
tribution of "host galaxy" redshift limits (275 objects) is 
shown as the blue hatched histogram in panel (c), also 
with an inset showing a blowup of z^ g < 0.7. The insets 
only go out to z — 0.7 because that is the maximum cali- 
brated redshift of the standard candle relation. For that 
reason, we also only compare the "host galaxy" redshifts 
to objects with z spec < 0.7 in the following. 

Even if one assumes that "host galaxy" redshifts are 
exact, the objects lacking spectroscopic redshifts are, 
on average, located at higher-redshift than objects with 
Zspec < 0.7. Considering the ~200 objects with reliable 
Zspec < 0.7, their average spectroscopic redshift is 0.33, 
compared to an average "host galaxy" redshift limit of 
0.48 for the 275 objects lacking spectroscopic redshifts. 
Later in this section we estimate that the "host galaxy" 
redshift limits are accurate to ±0.064. A Kolmogorov- 
Smirnov (K-S) test shows that the z spec distribution for 
objects with reliable z spec < 0.7 and the Zh g distribu- 
tion for objects lacking spectroscopic redshifts are statis- 
tically different; there is a negligible (10~ 24 ) chance they 
are drawn from the same parent distribution. There are 
likely some instances in our sample of extremely feature- 
less spectra that are nearby; in this case the AGN must 
be highly beamed, and the finite size of the SDSS opti- 
cal fiber limits the detectability of host galaxy spectral 
features. However, the majority of extremely featureless 
spectra are likely among the most distant, which could 
have implications for understanding BL Lac properties 
(for example, their cosmic evolution) if not accounted 
for properly. 

We check our "host galaxy" redshift algorithm by ex- 
amining 18 BL Lac candidates with lower limit spectro- 
scopic redshifts z spec < 0.7 (i.e., spectroscopic redshift 
limits derived from Mg II absorption). We do not see 
obvious host galaxy features in their spectra, so we ex- 
pect their "host galaxy" redshift limits to be smaller than 
(or approximately equal to) their spectroscopic redshift 
limits. This is indeed true: the average difference be- 
tween spectroscopic and "host galaxy" redshift limits is 
{z S pec - Zhg) = 0.18 (with z spec - z hg spanning -0.072 to 
0.41). These 18 objects are shown as red squares in Fig- 
ure [lOjs. While arrows designating lower limits are only 
drawn for the spectroscopic redshifts, the "host galaxy" 
redshift limits for these 18 objects should also be consid- 
ered lower limits. 

We also test our algorithm using a subset of 213 BL Lac 
candidates with reliable spectroscopic z spec < 0.7, and 
for which the host galaxy decomposition fit was pre- 
ferred over just a power law. Almost all of these objects 
have host galaxy to nuclear flux ratios larger than 0.2 at 
6581 A, so the above algorithm underestimates their host 
galaxy flux densities and should therefore overestimate 
their redshifts. These 213 objects are shown as black 
circles in Figure 110b . which confirms that the majority 
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2 3 4 5 
Redshift z spec 

Fig. 10. — Comparison between spectroscopic and host galaxy redshifts. This figure is in color in the online version, (a) The redshift 
distribution of all 723 BL Lac candidates (assuming redshift limits are exact). The red hatched histogram shows the redshift distribution 
of the 399 objects with spectroscopic redshifts, and the red solid histogram shows the 49 objects with lower limits from Mg II absorption. 
The inset shows a blowup of < z < 0.7. (b) Black circles: host galaxy redshift limits vs. spectroscopic redshifts for 213 objects with 
reliable z spec < 0.7. Host galaxy flux limits are estimated with Equation [4] which underestimates most host galaxies' true fluxes. Thus, 
z hg > Zspec for the vast majority. Red squares: 18 objects with lower limit z apec < 0.7. Host galaxies are not seen in their spectra, so the 
host galaxy redshift limits should be smaller than (or equal to) the spectroscopic redshift limits. The line z^ g = z apec is drawn to guide 
the eye. (c) Same as panel (a), except the blue hatched histogram shows the z^ g distribution for the 275 objects lacking spectroscopic 
redshifts. Objects lacking spectroscopic redshifts tend to be highcr-redshift on average, (d) Same as panel (b), except now the host galaxy 
flux estimates are exact (i.e., estimated with Equation [5j . The spectroscopic and host galaxy redshifts agree rather well (ctaz = ±0.064). 



(>90%) of the "host galaxy" redshift limit estimates are 
larger than their reliable spectroscopic redshifts. The 
line Zhg = z spec is drawn for guidance. 

Finally, a perhaps even better test is to estimate actual 
host galaxy flux densities for these 213 objects, and to 
then derive redshift estimates (instead of limits) . We es- 
timate the host galaxy flux densities (per unit frequency) 
at 6581 A for these 213 objects as: 



±0.064 (i.e., z > Zh g ± 0.064). We obtain similar uncer- 
tainties on the redshift limits if we add the uncertainties 
from our host galaxy flux measures to the rms error of 
the standard candle relation. We also note that the ac- 
curacy of our "host galaxy" redshift estimates is similar 
to uncertainties on "im aging" redshift limits placed by 
ISbarufatti et alj ([2005) for objects lacking host galaxy 
detections with HST. 



[fx,totai(Kff,R) - Aei«,o(1.067)- a *] , (5) 

where the second term inside the brackets uses the 
power law fit parameters (referenced at 6165 A) returned 
by Equation [3] to calculate the AGN's flux density at 
6581 A . We estimate the host galaxy flux density as the 
difference between the observed total flux density and 
the AGN flux density because it does not require a pri- 
ori knowledge of the redshift. We similarly do not use 
our knowledge of the spectroscopic redshift to estimate 
A{z). The spectroscopic and "host galaxy" redshifts 
match fairly well. They are compared in FigurefTOH. The 
line z spec = Z) lg is overplotted for guidance (i.e., the solid 
line is not a fit to the data). The "host galaxy" redshifts 
cluster along this line, with (Az) — (z spec — Zh g ) = 0.015 
and a dispersion of <r& z = ±0.064. 

In summary, we conclude from the above tests that our 
algorithm provides realistic redshift limits, and we esti- 
mate the "host galaxy" redshift limits to be accurate to 



5.2.2. AGN Component Optical Fluxes and Luminosities 

We also use the spectral fits to estimate optical fluxes 
and specific luminosities for the AGN component, ex- 
cluding the host galaxy flux when contamination is sig- 
nificant. For objects lacking spectroscopic redshifts, or 
for which the power law fit (Equation[2|) is preferred over 
the host galaxy decomposition (Equation [3]) , we assume 
the host galaxy does not contribute to the observed flux. 
In these cases, we estimate each AGN optical flux as the 
psf magnitude in the filter whose effective wavelength 24 
is closest to 5000 A rest-frame. These magnitudes are 
then corrected for G alactic extinction using the maps of 
ISchleeel et afl(fl998l) . We assume the "host galaxy" red- 
shift limit is exact to locate 5000 A rest-frame for objects 
lacking spectroscopic redshifts. 

If the host galaxy decomposition is preferred over the 
power law fit, then that source also has a spectroscopic 
redshift, which we use to identify the SDSS filter cover- 
ing (or closest to) 5000 A rest-frame. We use the AGN 
parameters from Equation [3] to estimate the AGN com- 

24 The SDSS filter effective wavelengths are 4686, 6165, 7481, 
and 8931 A for g,r,i, and z respectively. 
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ponent flux density at the chosen filter's effective wave- 
length. We then convert the decomposed AGN flux den- 
sities to optical AB apparent magnitudes, and we correct 
for extinction. The distribution of extinction corrected 
AGN component apparent magnitudes is shown in Fig- 
ure [9)3. Uncertainties on the decomposed AGN apparent 
magnitudes are typically around 0.03 mag. 

Next, we use the above apparent magnitudes to es- 
timate the AGN components' specific luminosities (per 
unit frequency) at 5000 A rest-frame (using each ob- 
ject's measured optical spectral index to perform the 
K-correction) . These luminosities are lower limits for 
objects with lower limit spectroscopic or "host galaxy" 
redshifts. 

5.3. Radio Matches in FIRST and NVSS 

We identify objects as radio emitters by correlating 
our 723 BL Lac candidates to the FIRST and NVSS ra- 
dio catalogs, and we find 561 matches within 2" to a 
FIRST source and 554 matches within 10" to an NVSS 
source. Our choice of a 10" search radius for NVSS is 
rather conservative, but we use it to avoid contamination 
from background radio sources (which increases rapidly 
for offsets r > 10".) There are 599 objects that match to 
either a FIRST or NVSS radio source. Among these 599 
objects, 516 objects match to both a FIRST and NVSS 
source, 45 objects match only to a FIRST source, and 38 
match only to an NVSS source. 

We expect the number of random mismatches between 
FIRST-SDSS and NVSS-SDSS to be small. Figure QT] 
shows cumulative distribution functions of separations 
between optical and radio source positions, Ar, nor- 
malized to the relevant search radii, r searc h = 2" and 
10" for FIRST and NVSS respectively. The distribu- 
tion functions are normalized to the total number of 
FIRST (561) and NVSS (554) radio matches. We ex- 
pect 90% of FIRST -SPSS matches to have Ar < 1" 
(|Becker et al.l I1995D . and we estimate that ^90% of 
NVSS-SPSS matches (for NVSS sources with flux den- 
sities larger than 5 mJy) should have Ar < 6" (see 
iCondon et al.1 fl998f ). The above 90% confidence radii 
(normalized to the relevant search radii) are overplotted 
as squares on each distribution function. Both distri- 
bution functions show that > 90% of FIRST-SPSS and 
NVSS-SPSS matches fall within Ar = 1" and 6" respec- 
tively. The horizontal dotted line in Figure [Tl] marks 
90% for guidance. 

We further investigate the chance of finding spurious 
matches to background radio sources by randomly offset- 
ting the SPSS source positions by ±3 — 10 times the rele- 
vant search radii for FIRST and NVSS. We then recorre- 
late the (offsetted) 723 BL Lac candidates to the FIRST 
and NVSS source catalogs. We repeat this test 10 times, 
and on average we only find 2 (1) spurious matches to 
FIRST (NVSS) radio sources (i.e., we expect < 1% of 
our sample to have spurious radio matches). Thus, our 
radio identifications appear to be statistically reliable. 

5.3.1. Radio Fluxes and Luminosities 

We take the radio flux density at 1.4 GHz to be the 
integrated flux density given in the FIRST catalog for 
all FIRST matches (561 objects). Because of FIRST'S 
better angular resolution, we use the FIRST flux density 




''search 



Fig. 11. — Cumulative distribution functions of the fraction of 
objects with projected source separations, Ar, between FIRST- 
SDSS (solid red line), NVSS-SDSS (green dashed-dotted line) and 
RASS-SDSS (blue dashed line) source positions. This figure is 
in color in the online version. The abscissa is normalized to the 
relevant search radii used for each survey, (2", 10", and 60" for 
FIRST, NVSS, and RASS respectively); the ordinate is normal- 
ized to the number of sources with matches in each survey (561, 
554, and 294 for FIRST, NVSS, and RASS respectively). The 
solid squares mark the separations we expect to include ~90% of 
matches from independent studies of source positional offset distri- 
butions. These distribution functions indicate that our radio and 
X-ray identifications are statistically reliable. 

even if the object is detected in NVSS. Since BL Lac ob- 
jects are compact, FIRST'S lower sensitivity to extended 
emission than NVSS (due to fewer short baselines) should 
not be important. If an object only matches to an NVSS 
source (38 objects), then we use its NVSS integrated flux 
density. If an object is not detected in either FIRST 
or NVSS but it is covered by the FIRST footprint (117 
objects), then we place radio flux density limits from 
FIRST as h.AGHz < 0.25 + (5cr rms ) mJy, where 
the noise in the FIRST survey at the object's sky coordi- 
nates (typically ^0.15 mJy beam -1 ). The 0.25 mJy term 
is a correctio n for the CLEAN bi as in the FIRST data re- 
ductions (see lBecker et al.l[l9 95). We do not place radio 
limits using NVSS (even for 7 objects outside the FIRST 
footprint) due to the increased potential for source con- 
fusion. The distribution of 1.4 GHz radio flux densities 
(and limits) is shown in Figure [9]:. 

We then calculate radio specific luminosities (per unit 
frequency) at rest-frame 5 GHz, assuming BL Lac spectra 
follow a power law with index a v = —0.27 (/„ ~ v~ av ) 
in the radio, which is the average radio spectral inde x 
of the 1 Jy BL Lac radio sample (Sti ckel et al.l fl99ll) . 
We again assume tentative and lower limit spectroscopic 
redshifts and "host galaxy" redshift limits are exact. 

5.4. X-ray Matches in RASS 

We identify 294 BL Lac candidates as X-ray emit- 
ters that match within 60" to a RASS so urce (in the 
RAS S Brig ht and Faint Source catalogs, iVoges et alJ 
119991 120001 ) . Given typical BL Lac SEPs and the flux 
sensitivities of RASS, SPSS, and FIRST, we expect most 
BL Lac candidates detected in X-rays by RASS to also 
be de tected in the radio (see Figure 1 of iPlotkin et al.l 
120081 ). All but 9 of the RASS detected BL Lac candi- 
dates are also detected in the radio by either FIRST or 
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NVSS. These 9 objects could be examples of extreme 
HBLs, they could have happened to be in low flux states 
during the FIRST and NVSS observations (or in high 
flux states during the RASS observations), or some of 
these 9 objects might be viable radio-quiet BL Lac can- 
didates, as even extant radio limits place about half of 
them in the radio-quiet regime (see m5.2\i . 

The cumulative distribution function of the RASS 
identifications is shown as the blue dashed line in Fig- 
ure HU We expect 90% of SPSS /RA SS source matches 
to have Ar < 25" (jVoges et al.1119991 ); this is indeed ob- 
served, as indicated by the blue square marking Ar = 
25" in Figure [TT] We also test our X-ray identifica- 
tions by randomly offsetting the SDSS source positions 
by ±3 — 10' and recorrelating the (offsetted) BL Lac can- 
didates to the RASS source catalog. After 10 iterations, 
wc find on average only 2 spurious matches (i.e., < 1% of 
our sample), further indicating that our 294 X-ray iden- 
tifications appear to be statistically reliable. 

5.4.1. X-ray Fluxes and Luminosities 

We estimate X-ray fluxes and specific luminosities (or 
upper limits) from the RASS count rates. We assume 
BL Lac objects have X-ray spectral indices of a v = 1.25 
(a value a pproximately intermediate between LB Ls and 
HBLs, see lSambruna et aT| [T996: Sambruna 1997), where 
the X-ray spectral index per unit frequency is again de- 
fined as f v ~ v~ a " . For the 429 objects lacking X-ray de- 
tections, we take the exposure time for the closest match 
in the RASS X-ray catalog. We then estimate upper 
limits for those objects' X-ray count rates as 6 photons 
divided by the ad opted exposure time. 25 We use the 
IStark et alj d 19921) hydrogen maps and the colden tool 
in CIAO (jFruscione et al.|[2006l ) to estimate the hydrogen 
column density along each object's sightline. We then 
use PIMMS to convert ROSAT PSPC X-ray count rates 
(or count rate limits) to unabsorbed broad-band (0.1- 
2.4 keV) fluxes (or limits). The distribution of broad- 
band X-ray fluxes and flux limits is shown in Figure [5Ji. 
We then convert the broad-band X-ray fluxes to X-ray 
specific luminosities (per unit frequency) at rest-frame 
1 keV, again assuming tentative and lower limit spec- 
troscopic redshifts and "host galaxy" redshift limits are 
exact. 

One X-ray undetected source in RASS (SDSS 
J232428.43+144324.3, z = 1.41) is weakly de- 
tected in our recent Chandra X-ray observation with 
0.0012 counts s" 1 from 0.5-6.0 keV. 26 We use its Chan- 
dra count rate to estimate it would have a RASS count 
rate of 0.0005 counts s~ x and an unabsorbed flux 1.25 x 
10~ 14 erg s -1 cm -2 from 0.1-2.4 keV. We assume a power 
law spectrum with photon index T = 2.25 and a c olumn 
density N H = 4.15 x 10 20 cm" 2 (|Stark et al.lll992T ) . The 
inferred RASS count rate is well below the limit for in- 
clusion in the RASS catalog. This source is outside of the 
FIRST footprint, and it is undetected by NVSS in the 
radio. It also remained undetected at 8.4 GHz in follow- 
on observations with the Very Large Array (VLA), with 

25 6 coun ts is the limit for inclusion in the RASS faint source 
catalog (sec Vogcs et al. 2000). 

26 The Chandra observation (Obsld 10386) was taken on 2009 
May 31 with the A dvanced CCD Imaging Spectrometer (ACIS; 
Garmirc et al. 2003) at the nominal S3 aimpoint. 



radio flux density /g.4 qhz < 0.3 mJy. 27 

5.5. Broad-band Spectral Indices 

Finally, we use the specific luminosities at each wave- 
band to calculate the broad-band spectral indices a ro 
and a ox (or limits in the cases of radio/X-ray non- 
detections). 28 We note that these values are relatively 
insensitive to the exact redshift used for their calcula- 
tion. The a ro — a ox plane is shown in Figure 1121 Ob- 
jects with RASS X-ray detections are shown as blue tri- 
angles, and objects lacking X-ray detections are shown 
as black circles. Red arrows designate a ro upper lim- 
its for objects lacking radio detections in FIRST/NVSS. 
All black circles are a ox lower limits, but for ease of dis- 
play we do not draw arrows for those limits in Figure 1121 
The blue square marks SDSS J2324+1443 (a ro < 0.003, 
a OI = 1.796) using its deeper Chandra X-ray and VLA 
radio observations. The dashed line marks the boundary 
between Low- and High-Energy Pe aked BL Lac objects 
(LBLs and HBLs, respectively, see iPadovani k. Giommil 
1995). 29 

The SDSS covers a large enough area to recover rare 
populations of objects, such as radio- weak AGN with fea- 
tureless spectra. The boundary between radio-loud and 
radio-quiet quasars is usually taken as a radio loudness 
para meter (i.e., radio to op t ical flux ratio) of R =10 
(see iKellermann et all Il989t iStocke et all I1992D . This 
approximately corresponds to a ro = 0.2 (marked with 
a dotted line in Figure [T2"]) , and we recover 86 feature- 
less spectra with a ro < 0.2 (not including the 7 objects 
outside the FIRST footprint.) We identify objects as 
radio-loud or radio-quiet assuming upper a ro limits are 
exact. Only 9 of the 86 radio-quiet objects have radio- 
detections, so the radio- weakness of many objects is even 
more extreme than we assume in this analysis. Sample 
spectra of 4 of these radio-quiet objects are shown in 
Figure [H 

Some of these radio-weak objects might be examples 
of beamed radio galaxies that simply populate the radio- 
faint tail of the much larger radio-loud SDSS BL Lac 
population. However, it is also possible that the SDSS 
is large enough to recover another rare population of ob- 
jects with featureless optical spectra. In this case, the 
radio-bright tail of the alternate population would over- 
lap with the radio-faint tail of the BL Lac population, 
making a secure BL Lac classification difficult with ex- 
tant data. The vast majority of these objects have spec- 
troscopic redshifts, and their average optical luminosity 
of vL v ~ 10 46 erg s" 1 (at rest-frame 5000 A; all 86 ob- 

27 This source was observed in the VLA D configuration on 
2007 March 30 (Program ID AP524). The VLA is operated by 
the National Radio Astronomy Observatory, a facility of the Na- 
tional Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 

28 The broad-band spectral index, for U2 > V\, is de- 
fined as ct u ,up = —log(Lv 2 /L vl )/log(u2/vi). Here, a ro = 

— log(L o /L r )/5.08, a ox = - log(L I /L o )/2.60, and a rx = 

— log(L I /L r )/7.68, where L r , L , and L x are the specific luminosi- 
ties (per unit frequency) at 5 GHz, 5000 A, and 1 keV respectively 
(Tananba um et aU ll979: Stock e~etaI1ll985l) . 

29 HBLs are commonly defined as a rx < 0.75 corresponding 
to F x /F r > IP - 11 with the X-ray flux in e rg s — 1 cm -2 and the 
radio flux in Jy (Padovani & Giommi 1995). The dividing line is 
also sometimes quoted as the similar F x /F r > 10 ~ 5 5 where both 
fluxes are in Jy (Wurtz 1994; Pcrlman et al. 1996). 
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Fig. 12. — Broad-band spectral indices a ro vs. ct ox . This figure is in color in the online version. The dotted line shows our adopted 
division between radio-loud (RL) and radio-quiet (RQ) AGN, and the dashed line shows the typical division (a rx = 0.75) between HBLs 
and LBLs. Red arrows show upper limits on a ro for objects lacking radio detections in FIRST/NVSS. Blue triangles: Objects that match 
to an X-ray source in RASS. Blue square: Object with a weak X-ray detection in a follow-on Chandra observation. Black circles: Objects 
that do not match to a RASS X-ray source. All black circles are actually lower limits on a ox , but to avoid cluttering the figure we do 
not illustrate these limits with arrows. There is a continuous population of BL Lac objects, with a smooth transition between HBLs and 
LBLs. However, the SDSS might be biased toward recovering excess HBLs. We also recover a large population of radio-quiet objects with 
featureless spectra; this population probably consists of a mix of BL Lac objects lying on the radio- weak tail of the parent distribution and 
another class of featureless objects like WLQs. Most of the recovered radio-quiet objects have z < 2.2, so they would constitute the largest 
sample to date of such exotic objects at these reds hifts. Some of the radio-quiet objects may populate a similar region of the a ro — ct ox 
plane as high-redshift WLQs (Shcmmcr ct al. 2009). 



jects have vL v > 10 43 erg s" 1 ) indicates that these are 
almost definitely AGN. We compare and contrast these 
odd objects with the larger radio-loud BL Lac sample in 

These 86 objects constitute an interesting population, 
whether or not they are not best described as BL Lac 
objects. For example, objects lying on the high-redshift 
tail (z > 2.2) might be classified as WLQs. Their emis- 
sion lines are plausibly intrinsically weak or absent rather 
than being washed out by s trong relativistically boosted 
iets dShemmer et all [20061 120091 iDiamond-Stanic et all 
2009). However, it is possible that some of these ob- 
jects still possess weak jets. Also, SDSS J2324+1443 (the 
blue square in Figure [T2]) suggests that deeper X-ray ob- 
servations of some of these objects may show them to 
have dissimilar X-ray properties compared to the radio- 
loud BL Lac candidates. We thus separate objects with 
a ro < 0.2 from the rest of the BL Lac candidates, and 
we refer to this population as weak-featured radio-quiet 
objects. 

5.6. Catalogs of Optically Selected BL Lac Candidates 

The final optically selected sample is presented in Ta- 
bles [5] - [HI Only the initial 5 entries of each table 
are shown here; the full versions are available in the 
electronic version of this article. Table [5] includes pri- 
marily empirical data from the SDSS database for the 



637 radio-loud BL Lac candidates. The IAU designa- 
tion, right ascension, and declination of each source are 
given in the first 3 columns. Our classification as a high- 
confidence ('H') or low-confidence ('L') BL Lac candidate 
is given in column (4). Column (5) indicates whether 
each source is identified as a BL Lac object in NED (as 
of May 2009): 'Y' means it is, and 'N' indicates it is 
not. Columns (6)-(10) give SDSS psf u,g,r,i,z band 
magnitudes. The spectroscopic redshift (z spec ) is given 
in column (11), and column (12) indicates if the red- 
shift is reliable ('R'), tentative ('T'), a lower limit from 
a Mg II absorption doublet ('L'), or unknown ('U'). The 
measured Ca II H/K depression, C, is given in column 
(13). Column (14) gives the AGN optical spectral in- 
dex a v measured in §5.21 Finally, column (15) notes 
whether a source is publ i shed a s an SD SS BL Lac can- 
didate in ICollinge et all (I2005L 'C 5'L lAnderson et all 
((20071 'A07') and /or IPlotkin et all (|2008|, 'P08'). This 
information is listed because not all SDSS BL Lac candi- 
dates are identified as such in NE D. We also note o bjects 
identified as z > 2.7 WLQs in [Shemmer et al.l (|2009L 
'S09') or in IDiamond-Stanic et all (I2009L 'D09'K TablelHl 
presents the same information for the 86 weak-featured 
radio-quiet objects. 

Table [7] includes primarily derived multiwavelength 
data for the 637 radio-loud BL Lac candidates. The 
SDSS source name and spectroscopic redshifts (z spec ) are 
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Fig. 13. — Sample spectra of 4 weak-featured radio-quiet objects with selected spectral features labeled. Absorption features are marked 
with dashed lines, and emission features are marked with dotted lines. No object matches to a R ASS X-ray source, and values for a ro 
are given in each panel. The spectral fits from the host galaxy decomposition described in £15.21 are overplotted as red lines. If there 
is significant host galaxy contamination, then the dotted red lines illustrate just the power law flux from the AGN. Top left: SDSS 
J103806. 71+382626. 5. This object has a tentative spectroscopic redshift assuming the strongest emission feature is Mg II. Although we 
assign this a tentative spectroscopic redshift, weak C II] may have also been detected. Top right: SDSS J110938. 50+373611. 6. This object 
has a reliable spectroscopic redshift and is detected in the radio by FIRST (with a ro = 0.168 placing it in the radio-quiet regime.) This 
may be an example of a BL Lac object lying on the radio-weak tail of the parent distribution. Bottom left: SDSS J132138. 86+010846. 2, an 
object with a reliable spectroscopic redshift. Bottom right: SDSS J233939. 48— 103539.3, a high-redshift object that shows the Lya forest, 
and weak N V and C IV emission. This object is typical of other SDSS high-rcdshift WLQs, and it is additionally listed in Table [3] as a 
new WLQ identification. 



TABLE 5 

Observed Parameters of Optically Selected Radio-loud BL Lac Candidates 



SDSS Name 


RA 


Dec 
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0.73824 


-0.41315 


H 


N 


21.10 


20.63 


19.89 


19.22 


18.85 


0.522 


R 


0.247 
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5.42608 
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H 


N 


20.15 
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0.648 


R 
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1.30 


C05;P08 


002200.95 + 000657.9 


5.50396 


0.11610 


H 


Y 
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TABLE 6 

Observed Parameters of Weak-featured Radio-quiet Objects 



SDSS Name 


RA 


Dec 


















Redshift 






(J2000) 


(.12000) 


(J2000) 


Classif. 


NED? 
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Flag 


C a u 


Comments 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


O) 


(10) 


(11) 


(12) 


(13) (14) 


(15) 


001741.84-105613.2 


4.42437 


-10.93703 


L 


N 


19.21 


18.98 
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R 
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18.89 
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T 
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025716.22-001433.5 
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L 


N 
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T 


0.61 




025743.73 + 011144.5 
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1.19572 


L 


N 
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18.76 
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1.712 


T 


1.01 




075331.84+270415.3 
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L 


N 
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18.90 
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0.74 
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repeated in columns (1) and (2) for easy comparison to 
Table [5] Column (3) includes our "host galaxy" redshift 
limits for objects lacking spectroscopic redshifts. An es- 
timate of the extinction corrected magnitude of just the 
flux from the AGN component is given in Column (4) 
(in the SDSS filter closest to 5000 A rest- frame). The 
logarithm of the corresponding AGN component specific 
luminosity at rest-frame 5000 A (in units of erg s _1 Hz -1 ) 
is given in Column (5). Column (6) includes a flag noting 
the radio survey from which we derive radio fluxes and 



luminosities. 'F' indicates the radio information comes 
from a detection in the FIRST survey; 'N' indicates the 
radio information comes from a detection in the NVSS 
survey; 'U' indicates that the source was undetected by 
both FIRST and NVSS, and our radio fluxes and lumi- 
nosities are upper limits derived from FIRST; 'X' means 
the object is outside of the FIRST footprint, and we do 
not include any radio information. Column (7) gives inte- 
grated radio flux densities at 1.4 GHz in mJy, and column 
(8) gives the corresponding logarithm of their specific lu- 
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minosities at rest-frame 5 GHz (in units erg s -1 Hz -1 ). 
The values in columns (7) and (8) are upper limits if 
the radio flag in column (6) is set to 'IP, and those en- 
tries are blank if the radio flag is set to 'X'. Column 

(9) contains a flag noting whether a source is detected 
('R') or undetected (TP) in the X-ray by RASS. Column 

(10) lists the broad-band unabsorbed X-ray flux from 
0.1-2.4 keV (in units of 10~ 13 erg s -1 cm -2 ), and the 
logarithm of the specific X-ray luminosity at 1 keV (in 
units of erg s _1 Hz -1 ) is given in column (11). Columns 
(10) and (11) are upper limits if the X-ray flag in column 
(9) is set to 'U'. Finally, broad-band spectral indices be- 
tween 5 GHz and 5000 A (a ro ), and 5000 A and 1 keV 
(a ox ) are given in columns (12) and (13) respectively. If 
the radio flag in column (6) is set to 'LP, then a ro is an 
upper limit, and a ox is a lower limit if the X-ray flag in 
column (9) is set to 'U'. The same information is pre- 
sented in Table [5] for the 86 weak-featured radio-quiet 
objects. 

There is some overlap between the phenomenological 
properties of BL Lac objects and WLQs. For example, 
6 of the serendipitous z > 2.2 WLQ discoveries listed in 
Table El additionally pass our BL Lac criteria (i.e., their 
emission features show measured REW < 5 A). Those 
6 objects are also included as weak-featured radio-quiet 
objects in Tables [5] and [S] 

6. DISCUSSION 

6.1. The Radio-Loud Subset 

The properties of the radio-lo ud subset are consis- 
tent with the results presented in iPlotkin et al.l (|2008f) . 
We refer the reader to that paper for a more detailed 
discussion: there is a smooth transition between LBLs 
and HBLs, and we also see a potential SDSS bias to- 
ward recovering HBLs. However, we note that many 
of the RASS X-ray flux limits are not sensitive enough 
to confirm this bias. Also, the radio-loud BL Lac can- 
didates tend to appear more luminous as the strength 
of the Ca II H/K break decreases (i.e., as one views 
the relativistic jet more directly along its axis, also see 
ILandt et al.l [20021. These properties are interpreted as 
generally supporting the standard orientation based uni- 
fication paradigm. 

6.2. What are the Weak-Featured Radio-Quiet Objects? 

The population of 86 radio-quiet objects with nearly 
featureless spectra is intriguing. Here we discuss the pos- 
sibility that these are a distinct population from the 637 
radio-loud BL Lac candidates. These objects are at sys- 
tematically higher redshift than their radio-loud coun- 
terparts. This is illustrated in Figure [T4"a. which shows 
the redshift distribution of all 723 BL Lac candidates. 
The dotted lines mark a region that is blown up in the 
inset, which shows the redshifts of the 86 radio-quiet ob- 
jects color-coded by their redshift. The blue histogram 
is for the 13 objects with z < 1. "Host galaxy" redshifts 
are shown as the black filled histogram for 5 radio-quiet 
objects lacking spectroscopic redshifts (all other redshifts 
for the radio-quiet objects are derived spectroscopically.) 
Mid-range redshifts (1 < z < 2.2) are shaded green, and 
large redshifts (z > 2.2) are shaded red. The open his- 
togram in the inset shows the distribution of the entire 
sample for reference. 



Only 5 radio-quiet objects lack spectroscopic redshifts, 
so there is not much concern for stellar contamination 
to the 86 object a ro < 0.2 subset. Plus, those 5 ob- 
jects lie outside of the gri box designed to remove stel- 
lar contaminants (see ij4.ip . as is shown in Figure [T4"b. 
Plotted are r — i vs. g — r colors, not corrected for ex- 
tinction or host-galaxy contamination; the color scheme 
is the same as in panel (a), with the 637 object radio- 
loud subset shown as black circles for comparison. The 
radio-quiet objects cover the same range of optical col- 
ors as the radio-loud BL Lac objects, further indicating 
significant overlap, at least phenomenologically, between 
the radio-loud and radio-quiet objects even if they are 
indeed two separate populations. It does appear, how- 
ever, that the optical colors of many of the 1 < z < 2.2 
radio-quiet objects tend to populate the bluer end of the 
radio- loud BL Lac distribution. But, their colors are typ- 
ical of normal spectro scopically confirmed SDSS quasars 
( Schneider et al. 20Q3) at similar redshifts that lack radio 
detections in FIRST (see contours in Figure [Mb) . 

We further explore the optical properties of the weak- 
featured radio-quiet objects by examining their optical 
luminosities as a function of redshift (Figure [14"b . which 
follows the same color scheme as panel b). The optical lu- 
minosities of the weak-featured radio-quiet objects follow 
the same trend with redshift as the radio-loud BL Lac ob- 
jects. Furthermore, the radio-quiet objects' luminosities 
are similar to the optical luminosities of ^65,000 normal 
spectroscopically confirmed SDSS quasars that lack ra- 
dio detections in FIRST. We thus conclude that, other 
than their lack of strong emission features, the weak- 
featured radio-quiet objects have similar optical proper- 
ties to both normal quasars and to radio-loud BL Lac 
objects. 

The radio-loud and radio-quiet objects also follow simi- 
lar optical flux distributions, as illustrated in Figure [l"4"h. 
The top panel shows the distribution of decomposed (ex- 
tinction corrected) AGN optical magnitudes for the en- 
tire 723 object sample, and a zoom in (with the blown 
up region marked by the dotted line) of the radio-quiet 
subset is shown in the bottom panel. There is a deficit 
of radio-quiet objects with mo, jn < 17, but overall the 
radio-quiet objects are not systematically fainter. The 
lack of very optically bright radio-quiet objects is likely 
due to a combination of their smaller sample size and 
their systematically higher redshifts. Thus, we can not 
attribute the lack of radio-detections for these objects to 
their having similar radio loudnesses as other BL Lac ob- 
jects but simply being optically fainter (and thus below 
the FIRST/NVSS detection limits.) In addition, even 
the radio limits placed by their non-detections firmly 
place them in the radio-quiet regime. 

If we assume the radio-quiet objects are "normal" 
BL Lac objects with comparatively weak radio emission, 
then we expect their radio emission to be weak because 
their SEDs peak in the UV/X-ray (i.e., most should be 
HBLs, or even extreme HBLs). We would therefore ex- 
pect a similar or even larger fraction of radio-quiet ob- 
jects to have X-ray detections in RASS compared to the 
radio-loud objects. However, around 40% of the radio- 
loud BL Lac candidates with mo :agn > 17 have X-ray 
detections in RASS, and only ~6% (5/86) of the radio- 
quiet objects have X-ray detections in RASS. Thus, al- 
though the recovered radio-quiet and radio-loud objects 
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TABLE 7 

Derived Parameters of Optically Selected Radio-loud BL Lac Candidates 



SDSS Namo 










Radio 


fl.4 GHz 




X-ray 










(J2000) 


spec 


*hg 


m 0,AGJV 


log £ 5000 A 


Flag 


(m.Iy) 


lo g L 5 GH 


Flag 




lo K L l keV 






(1) 




(3) 


(4) 


(5) 


(6) 


(7) 


(S) 


(9) 


(10) 


(11) 


(12) 


(13) 


000157.23-103117.3 


0.252 




18.17 


29.347 


F 


32.970 


31.822 


U 


3.811 


25.911 


0.487 


1.320 


000257.17-002447.3 


0.522 




18.80 


29.722 


F 


159.130 


33.146 


u 


3.373 


26.627 


0.674 


1.188 


001736.90+145101.9 


0.303 




17.73 


29.843 


N 


58.739 


32.236 


u 


6.946 


26.361 


0.471 


1.337 


002142.25-090044.4 


0.648 




18.46 


30.168 


F 


41.960 


32.752 


u 


3.976 


26.936 


0.509 


1.241 


002200.95 + 000657.9 


0.306 




19.21 


29.268 


F 


1.730 


30.714 


R 


21.140 


26.854 


0.285 


0.927 


NOTE. — Tablc^ is ava 


liable in 




rcty in the elc 


ctronic vcrsio 


n of this a 


rticlc. A porti 








rding its forn 


and cont 


ent. 



TABLE 8 

Derived Parameters of Weak-featured Radio-quiet Objects 



SDSS Name 
(J2000) 
(1) 


(2) 


z hg 1 
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n Q,AGN 
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log ^5000 A 
(5) 


Radio 
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(6) 


fl.4 GHz 
(m.Iy) 

(7) 


log L 5 GHz 
(S) 


X-ray 
Flag 
(9) 


F X log 
(10) 


L l keV 
(11) 


(12) 
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1.806 
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U 
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Fig. 14. — (a) The redshift distribution of all 723 nearly featureless spectra. The inset shows a blowup (marked by dotted lines in the 
larger panel) with z < 1, 1 < z < 2.2, and z > 2.2 radio-quiet objects colored blue, green, and red respectively. The black filled histogram 
shows "host galaxy" redshifts for 5 radio-quiet objects lacking spectroscopic redshifts. The radio-quiet objects are at systematically higher 
redshifts. (b) r — i vs. g — r, not corrected for extinction or host galaxy contamination; post- visual inspection color-cuts are overdrawn 
with dashed lines (see Q. Black circles: radio- loud BL Lac candidates; blue triangles, green crosses, and red star symbols: weak- featured 
radio-quiet objects with spectroscopic redshifts z < 1, 1 < z < 2.2, and z > 2.2, respectively. The optical colors of ^38,000 SDSS 
quasars with 1 < z < 2.2 that lack radio detections in FIRST are shown as contours (contours mark 100, 500, 1500, 2500, and 3500 
quasars). The weak-featured radio-quiet objects cover the full range of optical colors as the radio-loud objects, and the mid-redshift weak- 
featured radio-quiet objects in general have similar optical colors as normal 1 < z < 2.2 quasars. None of the radio-quiet objects lacking 
spectroscopic redshifts (shown as blue triangles circumscribed by squares) fall in the gri box aimed at removing DC white dwarfs (labeled 
WD), (c) Absolute magnitude vs. spectroscopic redshift. AGN component absolute magnitudes in the filter closest to 5000 A rest-frame are 
shown for the radio-loud and radio-quiet objects with spectroscopic redshifts (symbols are the same as panel b). Contours show ~65,000 
spectroscopically confirmed SDSS quasars without a match to a FIRST radio sou r ce {S chneider ct al. 2007); contours mark 100, 1000, 
2000, 3000, and 4000 quasars. Absolute magnitudes are taken from [Schneider ct al. (2007, i.e., in the i filter, assuming an optical spectral 
index per unit frequency of 0.5, and adopting the same cosmology as this study). The luminosities of the nearly featureless spectra are 
similar to normal quasars, indicating significant overlap in their optical properties (other than their weak or absent emission features), (d) 
Distribution of decomposed AGN optical apparent magnitudes. The top panel shows all 723 objects, and the bottom panel is a blowup 
(marked by the dotted line in the top panel). The distributions of radio-quiet objects at different redshifts are offset for clarity, and they 
follow the same coloring scheme as panel (a). 
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share similar optical properties, this hints that many of 
the radio-quiet objects might in fact be a distinct popula- 
tion f rom BL Lac objects (also see §4.1 of Shcmme r et all 
120091 who reach a similar conclusion for high-redshift 
WLQs). This, however, does not rule out that some of 
the radio-quiet objects indeed populate the radio-faint 
tail of the much larger radio-loud BL Lac distribution. 

We are unaware of any SDSS selection effects that 
would prevent us from recovering high-redshift radio- 
loud weak-featured objects or low-redshift radio-quiet 
weak-featured objects, if such types of objects actually 
exist. It is important to note, however, that some of the 
radio-loud BL Lac candidates lacking spectroscopic red- 
shifts have redshifts bounded by their "host galaxy" red- 
shift limits on the lower end, and z < 2.2 at the higher 
end. Thus, it is possible that many of the radio-loud 
BL Lac candidates actually do have redshifts similar to 
the bulk of the weak-featured radio-quiet objects. If that 
is the case, however, it is then surprising that we find so 
few radio-loud BL Lac candidates that show even weak- 
features at z > 1. While the redshift distribution of the 
radio-quiet population is not conclusive evidence on its 
own that we have actually recovered multiple popula- 
tions of weak-featured objects, it certainly is suggestive 
that these might constitute a distinct population from 
BL Lac objects (especially when taken in context with 
the unexpectedly small rate of X-ray detections for the 
radio-quiet objects.) The lack of a large parent popula- 
tion of radio-loud BL Lac objects at z > 2.2 does strongly 
argue against unifying the z > 2.2 weak-featured radio- 
quiet objects (i.e., W LQs) with BL Lac objects (e.g., see 
IShemmer et~aL|[2009l ). 

The weak-featured radio-quiet objects may populate a 
similar region of multiwavelength color space as high- 
redshift (z > 2.2) SDSS WLQs (see Figure 5 of 
IShemmer et al.l 120091 ) . Among the 86 weak- featured 
radio-quiet objects in Tables [S] and [51 13 have z > 2.2, 
of which 5 hav e prev i ously been identified as WL Qs 
(|Shemmer et al.l 120091: IDiamond-Stanic et all I2009D . 30 
We do not recover all previously known SDSS WLQs 
because of our stricter REW criterio n. The proposed 
z = 1.66 WLQ SDSS J0945+1009 IjHrvniewicz et all 
I2009D has Mg II emission with REW ~ 15 A, and even 
its emission is too large to appear in this sample. There 
are likely other similar objects in the SDSS database not 
included here. 

About two-thirds of the radio-quiet weak featured 
spectra are at 1 < z < 2.2. Have we discovered a large 
population of lower-redshift analogs to WLQs? Resolu- 
tion requires further observations to discriminate from 
radio-weak BL Lac objects, but either result is interest- 
ing. Polarimctric and variability monitoring, as well as 
searches for radio and/or X-ray emission would be useful. 
Even though the vast majority are not detected by RASS, 
their extant a ox limits do not yet exclude them from a 
regime that is typical for BL Lac objects. It would also 
be interesting to take ultraviolet spectra of the lower- 
redshift radio-quiet objects and infrared spectra of the 

30 Of the other eight z > 2.2 weak-featured radio-quiet objects, 
6 are also listed as new WLQ identifications in Table [3] another has 
a tentative redshift z = 2.22, and the remaining object has a lower 
limit redshift z > 2.21 derived from Mg II absorption. Neither of 
the latter two objects show an obvious Lyo forest, so they are not 
included as new WLQ identifications in Table [3] 



very high redshift radio-quiet objects to compare their 
emission at similar rest-frames. 

Regardless of their true nature, these radio-quiet ob- 
jects are a fascinating sub-population: such a large pop- 
ulation of radio-quiet BL Lac objects would pose a seri- 
ous challenge to the standard unification paradigm that 
maintains all BL Lac objects are beamed radio galaxies. 
Identification of these objects as lower-redshift analogs of 
high-redshift WLQs would greatly aid our understand- 
ing of those other exotic objects. Just the subset of 
radio-quiet featureless spectra with z < 2.2 would be 
almost as large as the number of z > 2.2 WLQs previ- 
ously identified by the SDSS (see IDiamond-Stanic et al.l 
2009). These complementary samples taken together 
should provide valuable clues to their true nature, as 
well as insight into the physical conditions necessary to 
form (or in this case not form) broad line regions near 
supermassive black holes. 

7. SUMMARY 

We present a sample of 723 objects with weak spectral 
features selected based on their optical properties. Ap- 
proximately 35% of the objects presented here are new 
discoveries, and ^75% were discovered by the SDSS. Be- 
low we summarize our main conclusions: 

1. Among the tens of thousands of contaminants re- 
moved during visual inspection, we note potentially in- 
teresting serendipitous discoveries, including 17 unusual 
white dwarfs, 9 unusual BALQSOs with extreme contin- 
uum dropoffs blueward of Mg II, and 20 high-redshift 
WLQs (with reliable z > 2.2 and Lya+N V REW < 
10 A). Six of these WLQs are additionally classified as 
BL Lac candidates because their emission features have 
REW < 5 A. 

2. Approximately 60% of our BL Lac candidates have 
spectroscopic redshifts (or limits from Mg II absorption). 
For objects lacking spectroscopic redshifts, we assume 
BL Lac host galaxies are standard candles, and we place 
redshift limits to those objects utilizing the fact that we 
do not detect host galaxy contamination to their SDSS 
spectra. We estimate their "host galaxy" redshift limits 
to be accurate to a z — ±0.064. We thus have redshift 
information (either from weak spectral features or limits 
placed from the lack of a host galaxy detection) for every 
single object in our sample. We infer the objects lacking 
spectroscopically-derived redshifts are at systematically 
higher redshift. We expect this result, since the objects 
with the weakest spectral features are arguably among 
the most highly beamed, which requires a low probability 
geometry. So, we are more likely to find the most highly 
beamed objects at higher redshift, where more volume is 
probed. However, there may also be some examples of 
highly beamed nearby objects in this sample as well. 

3. Each object emerges with homogeneous multiwave- 
length data coverage, and we include radio/X-ray fluxes 
(or limits) for all objects. Around 80% of the BL Lac 
candidates presented here match to a radio source in 
FIRST or NVSS, and we place upper limits on radio flux 
densities for the rest (except for 7 objects outside the 
FIRST footprint and not detected by NVSS.) Approxi- 
mately 40% match to an X-ray source in RASS, and we 
place upper X-ray flux limits for the rest. 

4. We use each object's high-quality SDSS spectrum 
to perform spectral decompositions of each object's host 
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galaxy and nucleus, and we report AGN component op- 
tical fluxes and luminosities. 

5. Based on radio fluxes and limits from the 
FIRST/NVSS radio surveys, we subdivide our sample 
into 637 radio-loud BL Lac candidates and 86 radio-quiet 
weak-featured spectra. There is significant overlap be- 
tween the optical properties of the radio-loud and radio- 
quiet objects. However, the radio-quiet objects tend to 
be at higher redshift, and they also have a smaller X-ray 
detection rate. 

We argue that the a ro < 0.2 objects are a mix of differ- 
ent populations: some of these objects (especially ones 
lacking spectroscopic redshifts or with small redshifts) 
probably populate the radio- weak tail of the much larger 
radio-loud BL Lac population. Other sources (especially 
the ones at higher redshifts) are likely WLQs. The very 
high-redshift nature of published SDSS WLQs is a selec- 
tion effect (because they are identified by the presence 
of the Lya forest), and this study provides a systematic 
search of the SDSS database capable of revealing lower- 
redshift analogs in large numbers. 

Even just the z < 2.2 radio-quiet objects constitute, to 
our knowledge, the largest sample of such weak-featured 
radio-quiet objects. The sample size is larger than many 
entire venerable BL Lac samples, and this subset is of 
similar size and complementary to the number of known 
high-redshift SDSS WLQs. 

6. The radio-loud objects in general support the stan- 
dard unification paradigm, which maintains that radio 
and X-ray selected BL Lac objects populate extreme ends 
of a single, larger, and continuous BL Lac population. 
We recover a mix of LBLs, IBLs, and XBLs, although 
the SDSS appears biased toward finding HBLs. How- 
ever, we note that the large observed scatter in BL Lac 
properties (e.g., the distribution of their AGN optical 
spectral indices) is probably due to innate differences as 
well. Even just the 637 radio-loud objects constitute one 
of the largest samples of BL Lac objects derived from a 
single selection technique. 

7. It will be interesting to assemble light curves for 
these BL Lac candidates once large-scale time-domain 
surveys, like the L arge Synoptic Survey Telescope (LSST, 
llvezic et al.ll20Q8t ) and the Panoramic Survey Telescope 
& Rapid Response System (Pan-STARRS 31 ), come on- 
line. It may even be possible to use flux variability from 
tho se surveys as a BL L ac selection criterion. For exam- 
ple, I Bauer et al.l (|2009fl recently found ^3000 objects in 
the Palomar-QUEST survey identified by their optical 

31 http: //pan- Starrs . if a. hawaii . edu 



flux variability; many of these objects may be blazars. 
Time-domain BL Lac searches over huge areas of the sky 
requires telescope access that is unrealistic with standard 
non-survey telescopes, but well within the capabilities of 
LSST and Pan-STARRS. 
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